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Abstract 



Over the past decades it has become evident that luminous matter amounts only to a 
small fraction of the energy density in the universe. More than 75 % is accounted for by 
what is called the dark energy, and about 20 % must exist in the form of some kind of non- 
relativistic dark matter. Among the candidates for the constituents of dark matter, the 
super symmetric (SUSY) neutralino is one of the most promising. This thesis discusses the 
search for supersymmetry at the future Large Hadron Collider (LHC) and the ongoing 
construction of one of the four large LHC experiments, the Compact Muon Solenoid 
(CMS), and focuses on the detection of signals from the annihilation of supersymmetric 
dark matter in the spectra of cosmic rays. 

CMS relies on the excellent performance of its components and thus requires strict 
quality control before their assembly. The final steps of assembly of 1061 silicon microstrip 
detector modules for the CMS tracker endcaps are performed at the 1. Physikalisches 
Institut B at the RWTH Aachen. A laser test facility for these modules was developed 
and is described in this thesis. In contrast to test procedures based only on the evaluation 
of pedestal and noise data, the test facility relies on the generation of signals in the silicon 
sensors by infrared laser illumination. Subsequent analysis of the signals allows reliable 
detection of module defects. The fully automatic test facility provides high throughput 
and easy operation for the series production of the modules. Its performance is validated 
by investigating a reference module with artificially prepared defects of three types: open 
wirebonds, short-circuited strips and pinholes. It is shown that all defects are clearly 
detected. In addition to defect detection, an indication for the type of defect is provided. 
In a further validation step, nine modules from a prototype series are investigated with 
the laser test facility. Confirming the earlier results on the reference module, defective 
strips are reliably identified. The results are in agreement with those from other test 
facilities using different techniques. 

Measurements of cosmic ray antiparticles, such as positrons, can impose strong con- 
straints on the nature of new physics beyond the Standard Model. However, cosmic 
ray positron measurements are experimentally very challenging due to the vast proton 
background. This thesis describes a novel approach of positron identification with the 
space-borne AMS-01 experiment, namely through the detection of bremsstrahlung con- 
version in a silicon microstrip detector. In contrast to earlier single-track analyses, this 
approach involves the selection and reconstruction of multi-track events. Subsequent to 
an introduction to cosmic ray physics and a description of the AMS-01 experiment, the 
discussion of the signal process shows that bremsstrahlung from protons is suppressed by 
a factor of more than 3-10^ with respect to positrons. The background to the positron 
sample can largely be suppressed using the topological and geometrical properties of the 
events. 

In order to obtain the highest positron selection efficiency possible, novel combinato- 
rial track finding algorithms were developed, particularly optimized for the signature of 
converted bremsstrahlung. By applying restrictions on the invariant mass of particles 
the background to the positron sample is largely eliminated. The remaining background 
contamination is determined from large samples of Monte Carlo data taking into account 
the effects of the geomagnetic field. It amounts to 26 % of the positron counts and is 



corrected for. In order to remove atmospheric secondaries from the positron and electron 
samples, a precise method involving trajectory backtracing in the magnetic field of the 
Earth was developed and is applied individually to all positron and electron candidates. 

The results of the positron measurement show that the bremsstrahlung approach ex- 
tends the sensitivity range of AMS-01 to positron momenta up to 50GeV/c, which is 
far beyond the original scope of the experiment. The precision of the positron measure- 
ment is statistically limited by the small number of particle counts. The positron fraction 

/{e^ + e~) is calculated for particle momenta in the range from 1 to 50GeV/c. For 
momenta up to 8 GeV it is found to be in good agreement with model predictions for 
background from purely secondary positron production, while at higher momenta there is 
indication for a positron overabundance. Therefore, the AMS-01 data lend further weight 
to the hints of a positron overabundance seen in the data from earlier experiments. 

In addition to the positron fraction, the absolute fluxes of positrons and electrons are 
calculated from the event samples of the present analysis. For this purpose, a method was 
developed which allows the determination of the geomagnetic transmission as a function 
of momentum and direction of incidence with high accuracy. The results of the flux 
calculation are found to be in very good agreement with earlier data, confirming the good 
performance of electron and positron selection with the bremsstrahlung approach. 

Finally, the positron fraction results from this analysis have been combined with re- 
sults from earlier experiments. In the combined results, the significance of the positron 
overabundance with respect to the background expectation for purely secondary positron 
production increases to 5.3 standard deviations, which would be reduced to 4.2 standard 
deviations without the results of the present analysis. Therefore, a statistical fluctua- 
tion causing the positron overabundance in the data with respect to the background-only 
expectation can now be excluded. 



Zusammenfassung 



Es ist eine der weitreichendsten Entdeckungen der letzten Jahrzehnte, daB die sichtbare 
Materie im Universum lediglich einen sehr kleinen Teil zu dessen Gesamtenergiedichte 
beitragt. Zu mehr als 75 % ist es von einer vollkommen unbekannten Form von Energie 
dominiert, fiir die der Name Dunkle Energie gepragt worden ist, wahrend etwa 20 % als 
gravitierende nichtrelativistische Dunkle Materie existieren muB. Unter den zahlreichen 
Kandidaten fiir die Konstituenten der Dunklen Materie ist das supersymmetrische (SUSY) 
Neutralino das vielversprechendste. Die vorliegende Arbeit befaBt sich mit der Suche nach 
supersymmetrischen Teilchen am zukiinftigen Large Hadron Collider (LHC) und dem 
Bau eines der vier groBen LHC-Experimente, dem Compact Muon Solenoid (CMS), und 
konzentriert sich auf den Nachweis von Signalen aus der Annihilation supersymmetrischer 
Dunkler Materie in den Spektren der Kosmischen Strahlung. 

Die einwandfreie Funktion seiner Komponenten ist eine wesentliche Voraussetzung fiir 
den Erfolg eines technisch so anspruchsvollen GroBexperimentes wie CMS und macht eine 
strenge Qualitatskontrolle wahrend ihrer Fertigung erforderlich. Die Endmontage von 
1061 Siliziumstreifendetektormodulen fiir den Spurdetektor von CMS wird im 1. Physi- 
kalischen Institut B der RWTH Aachen durchgefiihrt. Die vorliegende Arbeit beschreibt 
die Entwicklung eines vollstandig automatisierten Testverfahrens zur Qualitatskontrolle 
dieser Module. Es basiert auf mittels eines Infrarotlasers in den Siliziumsensoren an- 
geregten Signalen, deren Auswertung eine zuverlassige Erkennung von Defekten erlaubt. 
Um die Anwendbarkeit des Testverfahrens zu bestatigen, wurde ein Referenzmodul mit 
kiinstlich praparierten Defekten untersucht. Wie gezeigt wird, konnen alle Defekte ein- 
wandfrei detektiert werden, und es ergibt sich aus der Messung ein Hinweis auf den je- 
weiligen Defekttyp. Dariiberhinaus wurde das Testverfahren zur Untersuchung von neun 
Modulen einer Prototypserie verwendet. In Ubereinstimmung mit den Ergebnissen aus 
der Untersuchung des Referenzmoduls konnten Moduldefekte zuverlassig identifiziert wer- 
den. Die Ergebnisse werden von denjenigen anderer Testverfahren, die auf abweichenden 
Nachweistechniken beruhen, bestatigt. 

Die Energiespektren von Antiteilchen in der Kosmischen Strahlung, wie etwa von 
Positronen, konnen wichtige Hinweise auf die Natur der Dunklen Materie liefern. Die Iden- 
tifizierung von Positronen in der Kosmischen Strahlung stellt experimentell eine erhebliche 
Herausforderung dar aufgrund des groBen von Protonen verursachten Untergrundes. Die 
vorliegende Arbeit beschaftigt sich mit einem neuen Ansatz zur Identifikation von Positro- 
nen auf der Grundlage des Nachweises konvertierter Bremsstrahlungsphotonen in einem 
Siliziumstreifendetektor. Die Wahrscheinlichkeit der Emission von Bremsstrahlung von 
Protonen ist gegeniiber Positronen um einen Faktor von mehr als 3 ■ 10^ unterdriickt und 
erlaubt somit eine Unterscheidung beider Teilchensorten. Dieses Verfahren wurde anhand 
der Daten des weltraumgestiitzten AMS-Ol-Experimentes realisiert. Im Gegensatz zu 
friiheren auf Einzelspurmessungen basierenden Analysen erfordert es die Rekonstruktion 
von Mehrspurereignissen. 

Um eine bestmogliche Selektionseffizienz fiir Positronen zu errreichen, wurden speziell 
fiir die Signatur konvertierter Bremsstrahlungsphotonen optimierte Spurfindungsalgorith- 
men entwickelt. Schnitte auf die invariante Masse der gemessenen Teilchen eliminieren 
weitgehend den Untergrund zu den Positronkandidaten. Der verbleibende irreduzible 



Untergrund wurde mit Hilfe einer groBen Anzahl simulierter Ereignisse unter Beriicksich- 
tigung des Einflusses des Erdmagnetfeldes bestimmt. Die Zahl der Positronen wurde um 
diesen Untergrund, der 26 % der Signalereignisse entspricht, korrigiert. Zur Unterdriick- 
ung atmospharischer Sekundarteilchen unter den Positron- und Elektronkandidaten wurde 
eine Methode entwickelt, die auf der Riickverfolgung von Teilchenspuren im Erdmagnet- 
feld basiert. 

Die Ergebnisse dieser Positronmessung zeigen, daB der Bremsstrahlungsansatz den 
MeBbereich von AMS-01 bis bin zu Positronimpulsen von 50GeV/c erweitert, was weit 
jenseits der urspriinglich angestrebten Reichweite des Experimentes von 3,5GeV/c liegt. 
Die Genauigkeit der Positronmessung ist hierbei statistisch limitiert durch die Anzahl der 
gemessenen Positronen. Weiterhin wurde der Positronenanteil e"*"/ {e'^+e~) in der Kosmis- 
chen Strahlung fiir Teilchenimpulse im Bereich von 1 bis 50 GeV/c berechnet. Fiir Impulse 
bis bin zu 8 GeV/c ist dieser in guter Ubereinstimmung mit den Ergebnissen von Modell- 
rechnungen fiir den Untergrund von ausschlieBlich sekundar produzierten Positronen. Bei 
hoheren Impulsen hingegen zeichnet sich ein UberschuB von Positronen ab. Somit liefern 
auch die Daten des AMS-Ol-Experimentes Hinweise auf einen PositroneniiberschuB, wie 
er bei friiheren Experimenten bereits beobachtet wurde. 

Zusatzlich zum Positronenanteil wurden auch die absoluten Fliisse von Elektronen und 
Positronen aus der Anzahl der Kandidaten dieser Analyse errechnet. Hierzu wurde ein 
Verfahren entwickelt, das es erlaubt, die Durchlassigkeit des Erdmagnetfeldes fiir Elek- 
tronen und Positronen als Funktion des Impulses und der Flugrichtung der Teilchen mit 
hoher Genauigkeit zu ermitteln. Die Ergebnisse der FluBberechnung zeigen eine sehr gute 
Ubereinstimmung mit denjenigen friiherer Messungen. 

AbschlieBend wurde der im Rahmen dieser Arbeit ermittelte Positronenanteil mit den 
Ergebnissen anderer Experimente statistisch kombiniert. Aufgrund der Ergebnisse der 
vorliegenden Arbeit steigt die Signifikanz des Positroneniiberschusses gegeniiber dem er- 
warteten Untergrund von 4,2 auf 5,3 Standardabweichungen. Somit kann nun eine statis- 
tische Fluktuation als Ursache fiir den UberschuB von Positronen gegeniiber dem er- 
warteten Untergrund weitgehend ausgeschlossen werden. 
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1 Introduction 



The picture shown here is known as the Hub- 
ble Ultra Deep Field and was taken by the 
Hubble Space Telescope in 2004 p. It is the 
result of an ultra-long exposure - more than 
11 days effectively - to the light coming from 
a tiny area of the sky, as small as 3 by 3 arc- 
minutes. In other words, the image repre- 
sents far less than one millionth of the sphere 
that surrounds the Earth. An area of this size 
could easily be obscured by a grain of sand at 
the distance of an arm. 

Apart from a few faint stars which hap- 
pened to be in the field of view, about 10,000 
galaxies appear in the image. With the coarse 
but reasonable assumption that a typical galaxy 
contains about one hundred billion stars and that galaxies are uniformly distributed over 
large scales, one concludes that the total number of stars in the universe must be roughly 
of the order of 10^2, or 10,000,000,000,000,000,000,000. What an unimaginable amount 
of mass and energy fills the universe! 

Nevertheless, in the last few decades it became apparent that luminous matter and all 
the invisible dust and gas inbetween, whose elementary composition is in principle similar 
to that of matter in our familiar surroundings, amounts only to a surprisingly small 
fraction of the energy density in the universe. In reality, the universe must be regarded as 
essentially dominated by exotic and novel forms of matter, and an even bigger reservoir of 
strange non-material energy, whose nature yet remains a complete mystery. Consequently, 
the assertion of basically knowing what fills space must be abandoned and the terms dark 
matter and dark energy have been coined for the unknown bulk ingredients of the universe. 

At the same time the Standard Model of particle physics, which structures the knowl- 
edge about the elementary constituents and interactions of matter, has matured. It has 
been and still is extraordinarily successful in explaining practically all experimental obser- 
vations in particle physics and has received continuous confirmation with highest precision 
for more than two decades. Despite its success, it is widely assumed that the Standard 
Model must merely be a part of a superior and more comprehensive concept, i.e. it must 
be a low-energy effective theory truly valid only in the energy regime open to current 
experiments. Thus a large diversity of hypotheses has emerged to anticipate whatever 
may appear when the present experimental limits are exceeded. Among such hypotheses, 
the idea of supersymmetry is one of the most cogent and intriguing. 

It is becoming apparent that cosmology, astrophysics and particle physics are simulta- 
neously approaching the same problem from different points of view and that their synergy 
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1 Introduction 

may finally lead to a fundamentally new understanding of basic physics. In a manner of 
speaking, the universe acts as a huge autonomous reactor for all its constituents, even the 
possibly unknown, and offers by far the most powerful particle accelerators. Therefore, 
the study of cosmic ray particles can deliver strong constraints on those new phenomena 
which are as yet beyond the reach of any laboratory experiment and can point the way 
towards better theoretical and experimental approaches. With the upcoming collider fa- 
cilities, such as the LHC, it will thus be possible to target new discoveries with greater 
precision. 

This thesis discusses the search for supersymmetry at the future LHC collider and 
the ongoing construction of one of the four large LHC experiments, the Compact Muon 
Solenoid (CMS), and focuses on the detection of signals from super symmetric dark matter 
candidates in the spectra of cosmic rays. The second chapter introduces the question of 
the nature of dark matter together with an overview of the history and current state of 
the universe, as well as the evidence for the existence of as yet invisible new forms of 
matter. The supersymmetric extension of the Standard Model of particle physics leads to 
the presently most promising dark matter candidate, namely the hypothetical neutralino. 

The possible discovery of supersymmetry at the LHC is the subject of chapter three. It 
starts with an overview of the LHC facility and of the four large collider ring experiments. 
Subsequently, the discovery potential for supersymmetry of the LHC experiments is dis- 
cussed, including the production mechanisms and the decay channels of supersymmetric 
particles, as well as the experimental signatures which would allow their detection and 
the determination of supersymmetry parameters. The CMS experiment at the LHC is 
currently under construction and will feature the largest silicon microstrip tracker ever 
built. Such a complex detector relies on the excellent performance of its components and 
thus requires strict quality control before their assembly. Following an overview of the 
experiment and in particular of the tracker and its components, a laser test facility for 
silicon microstrip detector modules is described. The investigation of artificially created 
module defects establishes the applicability of the test facility and the results of tests of 
a prototype series of detector modules are presented. 

The fourth chapter addresses cosmic rays as possible messengers for signals from re- 
actions of supersymmetric dark matter particles. Subsequent to a description of the 
composition and characteristics of cosmic rays, their origin and propagation in the in- 
terstellar medium, as well as the modulation of cosmic ray spectra in the vicinity of the 
Earth are discussed. As it turns out, the annihilation of neutralinos from dark matter may 
leave its traces in cosmic rays in the form of an overabundance of particles, particularly of 
positrons, antiprotons, and gamma rays. The processes involved are therefore described 
and a characterization of the expected signatures in the corresponding particle spectra 
is given. The chapter ends with an overview of the evidence for neutralino annihilation 
which has so far been found in cosmic ray data. 

Such evidence still inconclusive and more experimental results are necessary to further 
investigate the hypothesis of supersymmetric dark matter. While new cosmic ray experi- 
ments are under construction or have recently started taking data, existing data from past 
experiments can be re-examined using new analysis techniques. The second part of this 
thesis reanalyzes the data taken by the space-borne AMS-01 detector in 1998 with the aim 
of determining the spectrum of cosmic ray positrons up to energies of 50 GeV, far beyond 
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the original scope of the experiment. After a detailed description of AMS-01 and its Space 
Shuttle flight in chapter five, chapter sixth introduces the challenge of cosmic ray positron 
measurements and the novel approach of positron identification through the detection of 
converted bremsstrahlung photons in silicon microstrip detectors. An in-depth descrip- 
tion of the analysis includes discussion of the signal and background processes, methods 
of particle track and vertex reconstruction, background suppression and correction for the 
background remaining among the signal candidates, as well as the complications which 
arise due to the circumstances under which AMS-01 was operated, namely the influence 
on particle detection of the spacecraft and of the magnetic field of the Earth. 

The results of the positron measurement are the subject of chapter seven, where the 
positron fraction - the ratio of positron to electron counts e~^/{e'^ + e~) in which several 
systematic uncertainties cancel - is presented as a function of particle momentum. It 
is shown that by using the conversion of bremsstrahlung photons to tag positrons and 
electrons in cosmic rays the AMS-01 data substantially improve the knowledge of the 
cosmic ray positron fraction up to the highest momenta so far accessible to experiments. 
To validate the tagging method, the chapter ends with the calculation of the absolute 
fluxes of positrons and electrons based on the results of this analysis and the comparison 
of these fluxes to the existing data [21 E]. 

The concluding chapter summarizes the results of this thesis. Finally, the positron 
fraction measurements from the most recent experiments [H E], El [7j are combined with 
the results of this analysis, giving the most precise positron fraction data yet available. 
The positron fraction measured in this analysis lends further weight to the hints seen in 
the data of [7] of a positron overabundance for momenta above 6 GeV. 
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2 Dark Matter in the Universe 



Over the past decades, numerous pioneering ideas and seminal discoveries have signifi- 
cantly broadened the understanding of the structure and evolution of the universe. The 
theoretical fundaments of space-time in an homogeneous and isotropic world based on 
General Relativity were already known in 1922 to contradict a static universe [Sj. Seven 
years later it was indeed discovered that the recessional velocity of distant galaxies is 
proportional to their distance, a fact which is nowadays referred to as Hubble's law [9]. 
Hence, the universe was found to be continuously expanding. The manifest idea of the 
Hot Big Bang - postulating that the universe primordially emerged from a tremendously 
dense and hot state - became increasingly popular, especially when it was shown that 
it could explain the observed abundances of the lightest nuclei fLOi [TT] . However, the 
strongest confirmation for the theory came with the discovery of the cosmic microwave 
background (CMB) in 1965 [12], an omnipresent electromagnetic radiation field with an 
almost perfect thermal 2.7 K black-body spectrum, which had been predicted before as a 
Hot Big Bang relic [13] . Recently, accurate measurements of the fluctuations in the CMB 
led to a precise determination of the cosmological parameters [H]. 

The first sky surveys and galaxy mappings surprisingly revealed that there is structure 
in the distribution of matter at very large scales in the order of 100 Mpc [I5]. Using the Big 
Bang model and assumptions about the types of matter that make up the universe, it is 
possible to predict the matter distribution. Furthermore, the survey of Type la supernovae 
yields a relation between the host galaxies' redshift and their distance, so that supernovae 
of this species are regarded as standard candles, allowing to investigate the expansion 
history of the universe [16] . The large-scale structure and supernova data as well as the 
CMB results clearly favor a class of models for which the term ACDM was coine cfl. The 
basic statement of the latter is that ordinary matter (also referred to as baryonic matter) 
constitutes only a small fraction of the energy density of the universe [T7] . More than 75 % 
is accounted for by what is called the dark energy, and about 20 % must exist in form of 
some kind of non-relativistic dark matter. The nature of these two most dominant forms 
of energy in the universe is one of the cardinal unsolved questions in modern cosmology, 
and the answer will doubtlessly give rise to fundamentally new physical concepts. While 
cosmology is presently far from a distinct approach to solving the dark energy problem, 
the nature of dark matter is possibly now close to being resolved. 

2.1 The History of the Universe 

The standard model of the universe is the Hot Big Bang model. According to this theory, 
the universe came into existence 13.7 billion years ago ^7\, starting in an enormously 
compressed and very hot state. Very little is known about the first 10~^^ s, the so called 

^CDM - cold dark matter; A denotes the cosmological constant. 
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2 Dark Matter in the Universe 




Fig. 2.1: A map of the universe: the distribution of distant galaxies exhibits structure on 
large scales [TH]. 



Planck Epoch, in which the three forces we know today (see § l2.3p were presumably unified. 
After the end of this period, the universe entered an approximately 10~^^ s long phase 
of exponentially quick expansion (inflation), during which its size increased by about 25 
orders of magnitude [I9] . The driving force of this process is believed to be a yet unknown 
scalar field, called the inHaton, whose initial quantum fluctuations from the Planck Epoch 
were blown up to macroscopic scales. These are nowadays represented by the large-scale 
structure (Figure I^TTl) . Simultaneously, a possible curvature of space-time in the early 
universe was stretched out by inflation, giving the universe its Euclidean geometry which 
is observed today [T7] . 

After the end of inflation, the energy in the universe was dominated by radiation and 
a sea of relativistic particle-antiparticle pairs near thermal equilibrium [20]. As the uni- 
verse further expanded and cooled, particle species dropped out of equilibrium once the 
thermal energy k^T fell below their mass, and vanished through annihilation. Accord- 
ing to Sakharov's conditions [21], the existence of CP- violating (see § 12.31) and baryon 
number violating processes out of thermal equilibrium may have lead to a tiny surplus of 
particles with respect to antiparticles which remained. At k^T ^ 200 MeV, 10~^ s after 
the inflation, a phase transition occured from a quark-gluon plasma to protons, neutrons, 
and pions, along with the leptons, antileptons, and photons (baryogenesis). When the 
universe was about 100 s old and k^T ^ 1 MeV, the temperature was low enough to allow 
for the creation of light nuclei such as D, ^He, ^He, and ^Li (nucleosynthesis). During this 
period, the last particle-antiparticle pairs, the electrons and positrons, annihilated. 

It took 370,000 years [Hj until the temperature was low enough {k-sT ^ 0.3 eV) for 
stable hydrogen atoms to be formed from the remaining protons and electrons (recombi- 
nation). At that time, the universe became transparent for radiation, which since then 
propagates almost freely. Today, it has cooled down to 2.7 K and is observable as the 
remarkably isotropic CMB. With the recombination and until the emergence of the first 
stars several hundred million years later, the Dark Age of the universe had begun. 
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2.2 Evidence for Dark Matter 




h"^ l<pc Radius (kpc) 

Fig. 2.2: The reconstructed surface mass density not associated with visible galaxies in the 
galaxy cluster CL0024 (white contours, in units of 2.1 • IO^Mqpc"^) superposed on an image 
by the Hubble Space Telescope [22] (left); rotation curve of the galaxy NGC6503 [53] with 
contributions from stars (disk), interstellar gas and the dark matter halo (right). 



The present universe can be characterized as spatially flat - Euclidean - and homoge- 
nous and isotropic on the largest scales. It is expanding so that objects move away from 
each other with a speed of 73km/s per Mpc distance [I7j. There is strong evidence that 
the expansion is presently accelerated [TB] , which can be explained by the presence of non- 
material dark energy with particular properties. More than 75 % of the energy density of 
the universe exists as dark energy, another 20 % as an unknown form of gravitating matter 
which otherwise does not observably participate in electromagnetic interactions. Energy 
densities are usually stated as fractions VL = p/pc oi the critical density pc which leads 
to a spatially flat universe. The value of pc corresponds to about six hydrogen atoms 
per cubic meter. In these terms, Q\ = 0.76 for dark energy, = 0.24 for the total 
matter density, and = 0.04 for baryonic matter only [T7]. The total energy density of 
the universe is = 1.003 ± 0.015 [H]. Matter is concentrated in galaxies, whose spatial 
distribution exhibits structure on scales of 1-100 Mpc in form of clusters, walls and voids. 
Since the hot and highly pressurized baryonic matter would have quickly washed out such 
structure, its formation can best be explained by the presence of dark matter. 



2.2 Evidence for Dark Matter 

Evidence for the presence of significant amounts of invisible matter is numerous and 
apparent from small to the largest scales. First indication was found as early as 1933 in 
the Coma galaxy cluster [21], which revealed a mass-to-light ratio two orders of magnitude 
larger than in the solar neighborhood. The mass of a cluster can be determined via several 
methods, including application of the virial theorem to the observed distribution of radial 
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2 Dark Matter in the Universe 




Fig. 2.3: Map of CMB temperature fluctuations by WMAP [5H]. 



velocities, by weak gravitational lensingj and by studying the profile of X-ray emission that 
traces the distribution of hot emitting gas [25]. Figure [2?2] (left) shows the reconstructed 
distribution of non- luminous - dark - excess matter in the galaxy cluster CL0024 [22], 
which apparently concentrates toward the cluster's center. Recent observations in this 
field are consistent with a contribution of 20-30 % from matter to the total energy density 
(see e.g. [2S])- This is way above what can be expected from the amount of visible matter 
alone. 

The most direct evidence for dark matter on galactic scales comes from observations on 
the rotation curves of galaxies. The radial velocity of a star at distance r from the center 
of the galaxy is (classically) given by: 



with G being the gravitational constant, and M(r) = / p{r)r'^dr the integrated matter 
density. Thus, the velocity should be falling oc 1/ ^/r outside the luminous galactic disk, as 
indicated by the dashed rotation curve in Figure [2^ (right) for the galaxy NGC6503 [23] . 
However, the data show that the actual radial velocity is largely a constant function of 
the distance from the center. This could be explained with the assumption that the 
galactic disk is embedded in a galactic halo of non-luminous dark matter, represented by 
the dash-dotted line in Figure (right). Basically all galaxies studied up to now share 
this behavior [27]. 

The spectacular recent CMB measurements from the Wilkinson Microwave Anisotropy 
Probe (WMAP) [11] and other experiments allow the determination of the cosmological 
parameters with unprecedented accuracy. The observed temperature anisotropies in the 



gravitational lens distorts the image of a background light source. In contrast to strong lensing, 
where there are easily visible distortions such as multiple images of the light source, weak lensing only 
leads to faint distortions on the percent level. 




(2.1) 
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Fig. 2.4 - A compilation of CMB power spectrum data and the best ACDM fit [281 
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Fig. 2.5 : Sensitivity of the CMB power spectrum to the baryonic and all-matter energy density 
contributions Ob (left) and (right) [29]. The parameter h equals about 0.74 [T7] . 



sky (see Figure [23]) are usually expanded in spherical harmonics Y^m (6', 0) as 

Sjrp +00 I 

^(^,0) = E E aemYem{e,<P) (2.2) 

1=2 m=-£ 

with the variance Ce = (|a^mP) [25]. Plotted as a function of £ in a power spectrum (see 
Figure [27^ . the variance reveals a noticeable structure with a series of distinct peaks. They 
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roughly correspond to resonances of acoustic oscillations in the photon-baryon plasma at 
the time of recombination. The oscillations occured around the equilibrium of radiation 
pressure and gravity. Since dark matter gravitates but does not interact with photons, 
it had a significant impact on the resonances. As a consequence, besides a large set of 
additional cosmological parameters, the CMB power spectrum is sensitive to the dark 
matter content of the universe. Figure 12.51 illustrates the effect of varying values of fib 
and fim on the power spectrum. The analysis of the data in the framework of the ACDM 
model [H] gives n^h'^ = 0.0223 ± 0.0008 and n^h^ = 0.126 ± 0.009 with h = 0.74 ± 0.03, 
which confirms that only about one fifth of the matter content of the universe is comprised 
of ordinary matter. 

As already stated, dark matter has left traces in the large-scale distribution of matter 
in the universe. Data from recent galaxy surveys [30, El] yield the spatial density distri- 
bution of large numbers of galaxies. The squared Fourier transform of this distribution 
- the power spectrum P{k) - constrains large-scale structure formation and cosmological 
parameters, such as the ratio fib/^m = 0.185 ±0.016 [HO], which is in excellent agreement 
with the CMB results. Furthermore, it can be inferred from the data that most of the 
dark matter particles must be nonrelativistic ("cold") and thus quite heavy, since rela- 
tivistic particles would have washed out structure on medium scales due to their large 
free-streaming length. 

Altogether, there is clear evidence for the existence of an unknown non-baryonic con- 
tribution to the matter content of the universe with an average density corresponding to 
roughly one hydrogen atom per cubic meter. Hence, the question for the nature of this 
dark matter arises. Is the current understanding of particle physics sufficient to explain 
its presence, or will new theoretical approaches be necessary? 



2.3 The Standard Model of Particle Physics 

The Standard Model of Particle Physics (SM) has, for many years, accounted for all ob- 
served particles and interactions (see e.g. [32] for reviews). In the SM, the fundamental 
constituents of matter are particles with half- integer spin (fermions): quarks and leptons. 
Their interactions are mediated by integer spin particles called gauge bosons. Strong inter- 
actions are mediated by gluons g, electroweak interactions by W^, and 7. The masses 
of the fermions and the and Z° are generated through the Higgs mechanism [33] via 
couplings to the Higgs boson H^. Gravitation is presently not covered by the SM, since 
General Relativity and the SM alone cannot be combined in one quantum theory. The 
left-handed leptons and quarks are arranged into three generations of SU{2)l doublets: 
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with the corresponding right-handed fields transforming as singlets under SU{2)l. Each 
generation contains two fiavors of quarks with baryon number B = 1/3 and lepton number 
L = and two leptons with B = and L = 1. Each particle also has a corresponding 
antiparticle with the same mass and opposite quantum numbers. The primed quarks are 
weak eigenstates related to mass eigenstates by the unitary Cabibbo-Kobayashi-Maskawa 
(CKM) matrix Vckm [3l[35] through 

( d'\ I Vud Vus Vub \ I d\ I d\ 

s' = V,d V,s W s \ = Vckm s (2.4) 

\b' J \Vu Vts Vtb I \ h ) \h j 

with complex parameters Vij. Under the assumption that exactly three quark generations 
exist and that the CKM matrix is unitary, the parameters can be constrained and partially 
absorbed into arbitrary phases. Four independent parameters remain, one of which can be 
interpreted as an imaginary phase. In the SM, this phase is responsible for the violation 
of the CP symmetry, which transforms particles into antiparticles with regard to their 
handedness. 

Local gauge symmetries play a fundamental role in particle physics. It is in fact in 
terms of symmetries and using the formalism of gauge theories that electroweak and 
strong interactions are described. The SM is based on the SU{?>)c ® SU{2)l ® U{1)y 
gauge theory, which undergoes the spontaneous breakdown 

SU{'i)c ® SU{2)l ® U{1)y SU{?>)c ® U{1)q , (2.5) 

where Y and Q denote the weak hypercharge and the electric charge generators, respec- 
tively, and SU{?>)c describes the strong (color) interaction, known as Quantum Chromo- 
Dynamics (QCD). This spontaneous symmetry breaking results in the generation of the 
massive and Z gauge bosons as well as a massive scalar Higgs field. 

The only viable dark matter candidate from the SM particle table is the neutrino v. 
The neutrino number density is known to be about 3/11 of the CMB photon density [HB] , 
thus yielding Vly = Y1, myC^ / (/;.^-93.8 eV) with the sum running over the three generations. 
With an upper limit on the sum of the neutrino masses of 0.72 eV/c^ (95% C.L.) from 
large-scale structure and CMB data [IT], the energy density contribution from neutrinos 
can be constrained to VLy < 0.014. Hence, neutrinos seem unlikely to constitute a rele- 
vant fraction of the dark matter in the universe. Further candidates from conventional 
physics include primordial black holes (see e.g. [36]) or massive object^ such as cold dwarf 
stars [37]. Up to now there has been no evidence that one of these is a good candidate, 
and in particular baryonic matter cannot constitute more than a small fraction of the total 
dark matter, thus it is likely that new physics needs to be introduced in order to approach 
the nature of dark matter. There is a wide range of candidates from non-Standard Model 
theories, including, amongst many others, axions and Kaluza-Klein states (see e.g. [211 ESj 
for reviews). However, the best motivated are provided by supersymmetric extensions to 
the SM. 



•^referred to as MACHOS - massive compact halo objects 
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Tab. 2.1 : Standard Model particles and their superpartners in the MSSM [25]. 



2.4 Supersymmetric Dark Matter 

Supersymmetry (SUSY) [38] is a generalization of the space-time symmetries of quantum 
field theory that transforms fermions into bosons and vice versa. In other words, the 
generators Q of supersymmetry act on particle fields as 

(5|fermion) = |boson) , Qlboson) = |fermion) , (2.6) 

such that there exist superpartners of the SM fields which mix or directly form supersym- 
metric particles (" sparticles" ) . Supersymmetry was originally postulated as a solution 
for the problem of quadratically divergent renormalization corrections to the Higgs boson 
mass, as such divergences cancel out with the introduction of equally heavy superpartners 
of the SM fields. Additionally, supersymmetry at TeV energies could modify the running 
coupling constants of electroweak and strong interactions in a way that they naturally con- 
verge to a common value at some very high energy scale ~ 2 ■ 10^^ GeV (GUtEI scale) [HU] . 
This physically attractive feature of SUSY is not present in the SM framework alone. 
Moreover, in contrast to the SM, the local extension of supersymmetry (supergravity) 
may provide a framework for the unification of the known forces with gravity. 

The Minimal Supersymmetric Standard Model (MSSM) has the smallest possible field 
content necessary to give rise to all the particles of the SM (see Table EH]). In the MSSM, 
unlike the SM, two Higgs doublets, corresponding to five physical Higgs bosons, are needed 
to give masses to all the fermions. The four neutralinos Xi 2 a 4 linear combinations 
of the superpartners of the gauge bosons B and and the neutral components of the 



"'GUT - grand unified theories 
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Higgs doublets (higgsinos), as 

X° = CuB + c,,W' + CsiH', + C4,^° , (2.7) 

while the two charginos xf,2 represent superpositions of the charged components and 
of the wino and the higgsinos. One additional ingredient usually used to constrain 
the MSSM is the conservation of R-parity. Ris a multiplicative quantum number defined 
as R = (^_l)3-B+i+25 ^Qj^ ^ particle with spin S, baryon number B and lepton number L. 
All SM particles have R = 1 and all their superpartners have i? = — 1. As a consequence 
of -R-parity conservation, SUSY particles can only decay into an odd number of SUSY 
particles, plus SM particles. Therefore, the lightest supersymmetric particle (LSP) is 
stable. 

If SUSY were an exact symmetry of nature, then particles and their superpartners 
would have identical mass. Since sparticles have not yet been observed, SUSY must be a 
broken symmetry. The simplest symmetry breaking mechanism, spontaneous symmetry 
breaking, is experimentally ruled out at a mass scale below 1 TeV |10]. Thus SUSY is 
thought to be broken spontaneously in a hidden sector, which consists of sparticles which 
do not directly interact with the visible SM particles or their superpartners. Symmetry 
breaking is then communicated to the visible sector via superheavy messenger particles. 
The cancellation of quadratic divergences mentioned above cannot be maintained in a 
broken SUSY, however the remaining divergences are only logarithmic and also propor- 
tional to the mass difference of particles and their superpartners. As a consequence, the 
sparticle masses are constrained and at least several of them should not significantly ex- 
ceed 1 TeV / c^ [3H] • Thus sparticles are in principle accessible to high energy colliders such 
as the LHC, which will be the subject of the next chapter. 

The MSSM has 124 free parameters, compared to 19 of the Sivlfl |10]. In the so called 
minimal supergravity (mSUGRA) models, gravitation is additionally taken into account, 
so that the spin-2 graviton has a spin-3/2 fermion superpartner called gravitino. Gravi- 
tational interactions then communicate the symmetry breaking from the hidden to the 
visible sector. In the mSUGRA case the number of parameters can be significantly re- 
duced by assuming that the model obeys a set of well motivated boundary conditions at 
the highest energies, thus eliminating all but 5 parameters: 

• tan /5, the ratio of the vacuum expectation values of the two Higgs doublets, 

• mi/2, the common gaugino mass at the GUT scale, 

• mo, the common scalar mass at the GUT scale, 

• Aq denoting the trilinear couplings of Higgs bosons to the sfermions, and 

• sign(/i), the sign of the Higgs boson mass parameter. 

This small number of parameters makes mSUGRA a fruitful field for supersymmetry 
phenomenology. Figure 12.61 (left ) shows the plane spanned by the parameters mo and 
mi/2 for fixed values of tan/3, Aq and sign(/i). Large regions can already be excluded 
on a theoretical basis, for example the region mo ^ fni/2, where the model does not 

^Hcre, the neutrinos are considered to be massless. 



13 



2 Dark Matter in the Universe 




Fig. 2.6: Constraints on the mSUGRA parameters niQ and mi^2 ^^^^ the neutrahno relic 
density in the allowed regions for tan/3 = 52, Aq = 0, > 0. The grey dotted and solid 
colored regions are excluded. The black contour denotes the cosmologically allowed region |41| 
(left); neutralino annihilation cross sections for parameters in the MSSM which are compatible 
with OcDM/i^ = 0.095 - 0.129 (dark points) or Ocdm^^ = 0.06 - 0.16 (light green points) for 
tan/3 = 50 02] (right). 



feature electroweak symmetry breaking (see previous section), or at low mo, where the 
LSP would be a charged sparticle. Further constraints are imposed by a wide range of 
data from e.g. WMAP, collider and cosmic ray experiments [H]. For the allowed regions, 
the corresponding neutralino relic density fl^h"^ is given. The black contour denotes the 
region consistent with results from WMAP and galaxy surveys. The allowed regions are 
typically located close to the borders of the excluded 

In mSUGRA and in a multitude of other SUSY models, the LSP is the lightest of the 
four neutralinos Xi- The mass of the x? is currently constrained by collider experiments to 
be larger than 46 GeV/c^ (95 % CL.) [32]. Neutralinos are Majorana particles, thus they 
can annihilate pairwise into SM particle states (see § l4.4p . The product of the annihilation 
cross section times the particle velocity aannV as a function of the Xi mass is displayed 
in Figure 12.61 (right) for those sets of MSSM parameters which are consistent with the 
WMAP observations with tan/5 = 50. Apparently, the data prefer neutralino masses of 
few hundred GeV and (Jg^nnV ~ 10~^^ cm'^ s~^. Given these numbers, it can be calculated 
that neutralinos must have left the primordial equilibrium at non-relativistic energies and 
that the cross section is indeed close to the rough expectation of ai^^^v ~ 10~^^ cm^ s~^ for 
particles in this mass range annihilating through weak interactions [13]. Taking all this 
into account, the neutralino appears as one of the most favorable candidates for particles 
which constitute large parts of cold dark matter in the universe: it has a high mass 
but does not carry charge, interacts only weakly with ordinary matter, has the correct 
abundance and non-relativistic speeds. Particles with these properties are 

generally referred to as weakly interacting massive particles (WlMPs). The fact that 
the neutralino does interact at all and, moreover, is embedded in a framework of other 
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sparticles, in principle opens the possibility to identify it as the constituent of dark matter 
by experiment. 

2.5 Approaches to Understanding Dark Matter 

It has been pointed out in this chapter that evidence for the existence of dark matter is 
widespread in astrophysics and cosmology. The fact that particle physics can provide re- 
liable models for dark matter candidates has established a new creative interplay of these 
domains in the recent past. A number of collider and non-accelerator experiments are 
currently operating or will soon be deployed that will shed light on dark matter. In addi- 
tion, there exist both direct and indirect non-accelerator dark matter search experiments 
that are ongoing or proposed [H]. Prospects for detecting dark matter and determining 
its properties are particularly bright in the case of the supersymmetric neutralino. 

From the view of particle physics, the dark matter question is strongly connected to 
those new physical phenomena which inevitably wait at energies beyond the reach of 
present experiments. Therefore, the particle physics program includes searching for su- 
persymmetric particles at new collider facilities, particularly the LHC, to possibly identify 
the LSP and determine its properties. The LHC and its experiments are currently under 
construction, so that their completion is naturally among the primary tasks regarding the 
search for dark matter. Since the field of supersymmetry contains an enormous variety of 
ideas, models and parameters, input from the domain of astronomy essentially helps to 
put constraints on supersymmetry and to focus the theoretical as well as the experimental 
work on realistic scenarios. 

From the view of astrophysics and cosmology, particle physics offers checkable pre- 
dictions about the behavior of dark matter and its relation to fundamental concepts of 
cosmology and the evolution of the universe. However, evidence for the existence of dark 
matter and insight in its properties come from various domains in astrophysics and cos- 
mology, and these can deliver data at energies which are currently unreachable for collider 
experiments. The so called direct detection of WIMPs is based on the measurement of 
nuclear recoils in elastic WIMP scattering processes [15], while the indirect search for 
dark matter is concerned with the detection of WIMPs through the analysis of their 
annihilation products in cosmic rays. At present, the experimental data show that the 
cosmic ray fluxes of positrons [7], antiprotons [IH), and gamma rays [17] are not quite in 
agreement with expectations, and that these discrepancies may well be interpreted as a 
signal of WIMP annihilation (see chapter HI). However, the data are still too sparse to 
make convincing statements. Hence, while waiting for new experiments to be deployed 
in the near future, putting more effort into the analysis of existing data can yield results 
with increased significance. 
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3.1 The LHC and its Experiments 

The Large Hadron Collider (LHC) [IHlHHllSO] is currently under construction at the CERN 
particle physics laboratory in the same underground tunnel where the LEP collider [5X1152] 
was previously housed. It will provide proton-proton collisions at a center-of-mass energy 
of 14TeV, as well as heavy ion collisions (lead-lead) up to 1.15 PeV. Two separated beam 
pipes will be used to circulate the protons or ions. The particles are forced on their 
trajectories by an 8.3 T magnetic field produced by superconducting magnets and are 
accelerated by superconducting radio frequency cavities. First operation of the LHC is 
scheduled for 2007 with beam energies of 450 GeV [53] . 

A certain process, which is characterized by its cross section a, occurs with a rate 
= Ca, where £. is the luminosity of the accelerator. The LHC is designed for a final 
luminosity C ~ 10^^ cm~^s~^ in the proton-proton mode, however it is planned to run at 
£ Ri 10^^ cm~^s~^ during the low luminosity phase for three years [M], which sums to 
approximately 30fb~^ integrated luminosity. To achieve the design luminosity the two 
beams will guide up to 2808 bunches of about 1.15 x 10^^ protons each. The bunch 
spacing will be 25 ns, which is equivalent to a bunch crossing rate of 40 MHz. Within 




Fig. 3.1 : The underground structures and experiments of the LHC [49| . 
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the total circumference of about 26.7km, there will be eight arcs and straight sections. 
Each straight section is approximately 528 m long and can serve as an experimental or 
utility insertion. As shown in Figure 13. four large experiments will be situated in the 
four beam crossing sites along the collider: 

• ATLAS [53], which is characterized by its sophisticated superconducting air-core 
toroid muon spectrometer, and features liquid argon calorimetry and silicon pixel 
and microstrip tracking, as well as transition radiation detectors, within a 2 T 
solenoid field; 

• CMS [56| [57], which is described in more detail below; 

• LHCb [SHI ESj, an experiment dedicated to studying phenomena related to CP- 
violation and rare B-meson decays at low luminosity, and 

• ALICE [60], which is customized for the investigation of ion-ion collisions. 

ATLAS and CMS are multi-purpose experiments and were optimized to perform high 
quality measurements in pp collisions of leptons, hadronic jets and high energy photons. 
The two detectors have been designed following similar guiding principles and differ mainly 
in the choice of detector technologies. While ATLAS has invested a large fraction of its 
resources into superconducting toroid magnets and into a set of precise muon chambers, 
CMS has put emphasis on the highest possible magnetic field combined with an inner 
tracker consisting solely of silicon pixel and silicon microstrip detectors, which provide 
high granularity at all radii. Main topics of the experiments' physics programs are the 
investigation of the mechanism responsible for electroweak symmetry breaking, which 
involves finding the Higgs boson, top quark physics and B-meson studies. However, since 
the LHC is often referred to as a discovery machine, a very important task is to search 
for new physics at high energies, such as SUSY. 

3.2 Discovery Potential of the LHC Experiments for 
Supersymmetry 

Assuming that supersymmetry is in fact realized in nature, its discovery at the LHC will 
be relatively straightforward [61]. There are many possibilities to create superpartners at 
a hadron collider: besides the quark-antiquark annihilation channel, there are numerous 
processes of gluon fusion, quark-antiquark and quark-gluon scattering. Gluon fusion leads 
to the largest cross sections, of the order of a few picobarn. SUSY would be revealed by 
an excess of events with a number of characteristic signatures over the standard model 
expectations. With i?-parity conserved, the final state in the decay chains of sparticles 
always contains LSPs, which escape the experiment undetected. Hence, signatures of 
SUSY are characterized by typically large values of missing transverse energy ^t- 

Table 13.11 gives an overview of the most common initial sparticle states and their ex- 
pected signatures at the LHC. At LHC energies, the total sparticle production cross 
section is dominated by strongly interacting gluinos and squarks. They initiate decay 
cascades which lead to the occurrence of final states with numerous jets and leptons and 
missing energy due to the existence of at least two LSPs and possibly neutrinos [61]. As a 
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Tab. 3.1: Compilation of initial sparticle states, main decay modes, and their signatures at 
the LHC [62]. 



consequence of LSP production, a complete mass reconstruction of gluinos and squarks is 
possible only in long decay chains. So SUSY signal observability is based on an excess of 
events of a given topology over known or expected backgrounds. Figure 13.21 (left ) shows 
the distribution of the quantity M^s - i-e. the sum of-^T and the transverse momenta pt 
of the four highest energetic jets - for events with > 100 GeV and at least four jets 
with Pt > 50 GeV/c, compared to the SM background. With an integrated luminosity of 
10fb~^, the signal over background ratio is approximately 10 for large values of Meg [S3]. 

The LHC can discover squarks and gluinos up to masses well in excess of 2 TeV/c^, which 
covers the entire parameter space over which supersymmetry can plausibly be relevant 
to electroweak symmetry breaking. In Figure [3^ (right) the reach of the gluino/squark 
search is displayed for 10 fb~^ of data. The solid lines show the 5a contours for final states 
with different lepton multiplicity. Additionally, contour lines for the resulting neutralino 
dark matter densities Qh"^ = 0.25, 0.5 and 1.0 are given. The cosmologically preferred 
density of Qh"^ ~ 0.1 (see chapter [2]) is well within the reach of the LHC. 

Sleptons can be produced in pairs, for example, directly via qq annihilation. If sleptons 
are more massive than the lightest charginos and neutralinos, they decay directly into 
these via — >■ ^^Xi2 ~^ ^ixt- Similar decay channels exist for sneutrinos. On the 

contrary, light sleptons can indirectly originate from charginos and neutralinos via e.g. 
X2 ~^ or xt ^ ^d^i where the charginos and neutralinos themselves dominantly 
stem from gluino and squark pairs. This leads to an experimental signature characterized 
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Fig. 3.2: Distribution of Mcs for the SUSY signal and the expected background for an 
integrated luminosity of 10fb~^ (left) [63]; reach {5a contours) of the gluino/squark search in 
various lepton multiplicity final states (OS - opposite charge sign; SS - same sign) for 10fb~^. 
The neutralino dark matter density contours of Qh'^ = 0.25, 0.5 and 1.0 are also shown. Dash- 
dotted lines denote squark/gluino isomass contours (right) [61]. 



by 2 isolated leptons, missing transverse energy, and the absence of jets, which allows to 
detect sleptons up to masses of about 350GeV/c^ pTj . 

Following the possible detection of SUSY signals in channels such as those described 
above, the most important tasks would be to measure the sparticle masses and to narrow 
down the range of model parameters. Leptonic decays of the second-lightest neutralino, 
X2, have a useful kinematical feature: the two- lepton invariant mass spectrum has an edge 
near the kinematical upper limit with a maximum value of m'^l- = m^o — m^o in the 
case of direct three-body decays via X2 ~^ and 

mfn- = \j ("^|o - mf) (mf - m|o)/m,- (3.1) 

in the case of two-body cascade decays X2 ~^ i'^i^ i^i'xl, thus allowing a mea- 
surement of the mass difference of the two neutralinos. The X2 dominantly stem from 
gluinos, which decay into heavy quark-squark pairs, as for example g bb, b "x^b. 
Thus the experimental signature is characterized by the occurrence of isolated same-flavor 
opposite-charge lepton pairs in conjunction with at least four jets. Figure (left) shows 
the invariant mass of lepton pairs from events which match these requirements for an 
integrated luminosity of 30 fb~ and a particular choice of parameters. The background 
expected from SM processes is also shown. While the narrow peak around 90GeV/c^ 
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Fig. 3.3 : Distribution of the dilepton invariant mass for mo = 800 GeV/c^ ,"^1/2 = 200 GeV/c^ , 
Aq = 0, tan/3 = 10, /U > 0, and the SM background, for vCint = 30fb~^ f^eftj; same with mixed- 
flavor pair distribution subtracted (right). Taken from [55] . 



stems from Z° production in the decays of heavy neutralinos, the edge at lower invariant 
masses indicates the mass difference Am^ = m^o — m^o , which in this case approximately 
equals 68 GeV/c^. The background is mostly due to it or VTiy+jet production. However, 
in that case the final state contains as many same-flavor leptons as different-flavor ones 
and with identical distributions. Hence, by subtracting the different-flavor distribution, 
the SM background can be canceled up to statistical fluctuations (Figure [3T3l right). Am^ 
can then be determined with a statistical accuracy of about 50MeV/c^ [53] . 

As eq. (13.11) suggests, the dilepton edge is also sensitive to the intermediate slepton 
mass. Through analysis of particular kinematical distributions, the masses of the two 
lightest neutralinos as well as the slepton mass can be determined [61]. Additionally, 
events near the edge can be used to reconstruct the X2 momentum vector. This would 
permit the search for resonance structures in the distribution of the X2+j^^ invariant mass 
and possibly provide access to the mass of the decaying squarks. Apart from these event 
topologies, a variety of other signatures would facilitate the determination of sparticle 
masses and SUSY model parameters [55], |57], [61], (63] [HH [65|. The ultimate goal of such 
studies would be to use very many measurements to make an overconstrained fit to the 
model, rather in the same way that current data are used to test the SM. 



3.3 Silicon Microstrip Detectors for CMS 

3.3.1 The CMS Experiment 

CMS is a general purpose high energy physics experiment which is particularly optimized 
for the detection of particles stemming from proton-proton collisions. Its design goals 
comprise, amongst others, providing a high-quality and redundant muon system, best 
possible electromagnetic calorimetry and excellent tracking in a strong magnetic field. 
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Its layout is that of a classical hermetic detector with a barrel region and two endcaps, 
covering practically the full An solid angle. CMS will have a length of 21.6 m, a diameter 
of 14.6 m, and 14500 1 of weight [66]. Figure [331 gives an overview of the experiment and 
its subdetectors, which will briefly be discussed in the following. 

The Magnet System [67] consists of a superconducting solenoid coil in conjunction 
with a saturated iron yoke for flux return. The 12.5 m long solenoid will provide an 
homogeneous magnetic field of 4 T in the central region, which decreases to about 1.7 T 
in the outermost parts of the experiment. 

The Muon System [68] is integrated in the return yoke of the magnet system and will 
be used to identify muons and measure their charge sign and momentum. Furthermore, 
it will play an important role for triggering the CMS experiment. The Muon System 
consists of three subdetectors: drift tube chambers in the barrel region, cathode strip 
chambers in the endcaps, and fast resistive plate chambers in both the barrel and end- 
cap. When combined with data from the tracker, the expected momentum resolution 
^pt/Vt foi' muons ranges from 1 % for < 10 GeV/c to about 6-17% for = 1 TeV/c, 
depending on the direction of the emerging particle. The Muon System also provides a 
track reconstruction efficiency of better than 90 % in the transverse momentum range up 
to 100 GeV/c, decreasing to > 70 % at = 1 TeV/c. 

The Calorimeters [691 170] of CMS will measure the energy and direction of elec- 
trons, photons and jets and, moreover, deliver information crucial for the trigger system. 
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Fig. 3.5: Cross section through one octant of the CMS siUcon strip tracker |72] . The pixel 
detector in not shown. 



They are divided into two parts, the electromagnetic calorimeter (ECAL) and the hadron 
calorimeter (HCAL). The complete ECAL and most parts of the HCAL are located within 
the magnet coil. The ECAL is composed of 82728 scintillating PbW04 crystals with a 
depth of 25.8 (barrel) and 24.7 (endcap) radiation lengths, read out by avalanche photodi- 
odes and vacuum phototriodes, respectively. In the barrel part, the ECAL granularity is 
2.2 X 2.2 cm^, providing an excellent separation power for close-to-collinear particles. The 
endcaps have less fine granularity and a silicon strip preshower detector will be installed 
in front of them to help with 7r°-7 separation. 

The HCAL consists of a barrel and two endcap (HC) sections inside the magnet coil 
and a tail catcher outside the magnet. It is a sampling calorimeter composed of copper 
absorber plates interleaved with 4 mm (front part) or 8 mm (back part) plastic scintillator 
tiles read out by hybrid photodiodes. In addition two forward quartz fiber calorimeters 
(HP) are placed around the beam pipe beyond the endcap magnet yokes to provide her- 
metic closure of the detector. The energy resolution as/E of the calorimeters is given 
by 

with E in units of GeV and a = 2.7%, a„ = 155 MeV, b = 0.5% for the ECAL (barrel, 
T] = 0), and a = 100%, b = 4.5% with negligible ct„ for the HCAL with the ECAL in 
front. 

The Tracking System [711 173] consists of two subsystems, the pixel detector and 
the silicon strip detector, with an overall transverse momentum resolution of 6pt/pt ^ 
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Fig. 3.6: Schematic view of the fourth TEC disk (z ~ 170 cm). It carries 16 petals - eight 
on the front and eight on the back side - on which the sihcon strip modules are mounted. The 
petals on this disk do not carry modules of ring 1 (compare Figure 1331) . 



(15 -px/ TeV©0.5) %. The pixel detector forms the innermost part of the tracker, close to 
the beam pipe. It consists of three 53 cm long barrel layers and two endcap disks on each 
side of the barrel. The barrel layers are placed at distances of 4.4, 7.3 and 10.2 cm from the 
beam line, while the two disks are situated at 34.5 and 46.5 cm from the interaction point. 
To achieve an equally good hit resolution of 15 /im in the transverse and the longitudinal 
planes, a design with a square pixel shape of dimensions 150 x 150 /xm^ and thickness 
300 yum is used 

The pixel detector is enclosed by the silicon strip detector (SST). With more than 15000 
silicon modules, adding up to a total active silicon area of about 200 m^, and roughly 10 
million electronics channels, it will be the largest silicon detector ever built [72]. Figure 
displays a cross section through one octant of the SST and shows its division into four 
subsystems: the Tracker Inner Barrel (TIB), the Tracker Inner Disks (TID), the Tracker 
Outer Barrel (TOB) and the Tracker Endcaps (TEC). The TEC and their components 
are described in the next section. 

3.3.2 Silicon Microstrip Detectors 

Each TEC consists of nine circular disks perpendicular to the beam direction, on which 
the the silicon sensors are mounted via modular support structures (petals) as shown 
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Fig. 3.7: Photograph of a fully assembled TEC petal (b-side), carrying modules of rings 2, 4 
and 6, and a view of a 22.5 cm long ring 6 module |72j . 



in Figure 13.61 One or two silicon sensors are integrated in modules together with the 
frontend readout electronics. The modules are arranged on both sides of the support 
structures forming up to seven concentric rings around the beam pipe with the sensitive 
strips oriented radially. Modules belonging to odd-numbered rings (1, 3, 5 and 7) are 
mounted on one petal side, while those of even-numbered rings (2, 4 and 6) are mounted 
on the other side, so that there is an overlap between adjacent rings in r and 0. Rings 3, 
4, 6 and 7 have single sided modules consisting of one single silicon plane. Modules of the 
rings 1, 2 and 5 are double sided (stereo) with two coplanar silicon planes mounted back 
to back and tilted at an angle of 100 mrad with respect to each other. Since at a given z 
position the silicon microstrip sensors are sensitive only to one of the coordinates spanning 
the sensor plane, single sided modules yield measurements of only the (p and z coordinates 
of particle hits, while double sided modules measure all coordinates simultaneously. The 
spatial resolution of single hit measurements in a TEC sensor is about 30 fim. Figure 13.71 
shows a fully assembled petal (b-side), carrying modules of rings 2, 4 and 6. Due to the 
overall circular geometry of the TEC the modules have a trapezoidal shape and increase 
in size towards the outer edge of the disks. 

For the modules of ring 6, the Physikalisches Institut IB at the RWTH Aachen is 
involved in the final assembly steps of wirebonding. Such a ring 6 module is shown 
in the enlarged view on the right side of Figure 13.71 Two carbon fiber frame legs are 
connected via a graphite plate (cross piece), serving as the mechanical support structure 
for the two trapezoidal silicon sensors named W6A and W6B. The latter are electrically 
insulated from the support structure with bands of kapton foil, which at the same time 
carry electrical circuitry for the application of high voltage to the sensors. The frontend 
electronics, assembled on a kapton/copper board laminated on thin ceramics (hybrid), 
are glued on the cross piece. A thin glass plate with strip conductors, the pitch adapter, 
guides the signals from the sensors to the readout electronics. Each of the single sided 
(double sided) sensors has 512 (768) readout strips. For ring 6 the strip pitch varies from 
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163 to 205 /iin. The inter-strip connections from sensor to sensor, sensor to pitch adapter, 
and pitch adapter to hybrid are estabhshed via 25 /im aluminum wirebonds. 

The structure and working principle of a silicon microstrip sensor is illustrated in Fig- 
ure I3.8[ Strips of p+ doped silicon are implanted into the surface of the 500 /im thiclfl n^ 
doped silicon bulk. They are insulated via a continuous Si02/Si3N4 layer, and are cov- 
ered with aluminum metallization for readout purposes (strips). The strip width varies 
between 53 and 67 fim for ring 6 modules. In this configuration, the implants form a series 
of p-n junctions, with the aluminum strips capacitively coupled to them. On both ends 
of each strip, the aluminum broadens and forms two AC pads to house the wirebond. In 
addition to the AC pads, DC pads are located beyond the ends of each strip. These pro- 
vide a direct connection to the implants and are thus not used for wirebonding. The DC 
pads are connected to the bias ring via 1.5MQ polysilicon bias resistors, hence providing 
all implants with a common potential. The opposite surface of the sensors is finished with 
a thin n"*" silicon layer and finally coated with an aluminum backplane. 

The frontend electronics on the hybrid are equipped with APy25-Si§ [TU] integrated 
circuits (IC). The APV25-S1 is a 128 channel analog pipeline with 192 columns of analog 
storage. Each strip is connected to one single electronics channel. Strip signals are 
amplified into 50 ns shaped pulses of magnitude 100 mV per 25000 electrons. These are 
sampled at a rate of 40 MHz and stored in the pipeline. Useful data are marked after a 
programmable latency, and held in the pipeline until such a time that they can be read 
out. If no trigger is received within 4.8 fis the pipeline cells are overwritten. A set of 512 
strip signals is referred to as an event in the following. The APV25-S1 can be operated 
in two different readout modes; in peak mode only one sample per channel is read from 
the pipeline while in deconvolution mode three samples are read sequentially, and the 
final output is the weighted sum of all three. In the latter case, the shaping constant is 
reduced from 50 ns to 25 ns at the expense of slightly higher noise. Besides a set of other 
parameters, the gain of the amplifier stage of the APV25-S1 can be varied. In addition 

^For rings 5-7; sensors of the rings 1-4 have a thickness of 300 /xm instead. 
^APV - Analogue Pipchnc (Voltage mode) 
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to the APV25-S1, the hybrid houses an IC to multiplex a pair of two APV25-S1 analog 
outputs onto a single differential line (MUX), the detector control unit (DCU) IC for slow 
control purposes and an IC for processing the trigger signal supplied by the CMS trigger 
electronics (TPLL) [77]. 

With a voltage V applied to the sensor backplane and the bias ring grounded, the p-n 
junctions are reverse biased, and depletion zones form at their locations. The width d of 
the depletion zones can be calculated in good approximation as 



where A'd is the effective doping concentration and q is the elementary charge. Full 
depletion is reached if d equals the sensor thickness ds, hence the corresponding depletion 
voltage Vdep is found to be 

Vaep«^|. (3.4) 

equal to about 100 V for A'^d ~ 0.5 ■ 10^^/cm~^ [7T] in a 500 /im thick ring 6 sensor. 
However, it must be pointed out that Nd and consequently Vdep will strongly vary as the 
silicon is irradiated during operation of the CMS experiment, essentially leading to a type 
inversion from the initially n-type silicon to p-type silicon |78j . 

Charged particles traversing the sensor ionize the silicon and create electron-hole pairs. 
A minimum-ionizing particle (MIP) deposits a most probable energy of 260 eV within 
1 /im of material along its flight path [32]; with an energy of 3.6 eV required to create 
one electron-hole pair in silicon [7^, it follows that a MIP releases about 36000 pairs or 
5.8 fC per charge sign when traversing a ring 6 sensor. The charge carriers are separated 
in the electric field of the depletion zones and move towards the electrodes. The signal is 
then created on the aluminum strips by induction and carried to the frontend electronics 
through the wirebonds. 

Irradiation with light also generates electron-hole pairs in silicon. The penetration 
depth a - i.e. the depth at which the intensity / of the incident light has decreased to 
1/e of its initial value Jq - strongly depends on the wavelength A. In the near infrared 
domain, it is a monotonously increasing function of A, and reaches a value of about 400 fim 
at A = 1050 nm in the case of pure silicon [SU] • The mean depth x of energy loss of light 
in a silicon sensor of thickness b as the expectation value of the intensity distribution 
I{x) = Iq ■ exp(— x/a) along x is 

X = a — b ■ rr ■ (3.5) 

For A = 1050 nm and b = 500 /xm, x equals 200 /im, which is close to x = 250 /xm obtained 
for the uniform energy loss attributed to a MIP. Although the silicon of the sensors is not 
pure but doped, the effect of the doping concentration on x is negligible for photons with 
energies above the 1.12 eV band gap of silicon, equivalent to a wavelength of 1107 nm. 

There exist several sources of noise superimposed on the strip signals. Besides noise 
due to thermal excitations in the silicon, to the readout electronics and low-frequency 
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voltage noise, the module output suffers to a small extent from the so called common 
mode noise (CMN), which is basically caused by electromagnetic pick-up at the APV25- 
Sl preamplifier inputs. One characteristic feature of the CMN is that it simultaneously 
affects groups of strips, and thus appears as a shift of their common baseline randomly 
varying with time. The APV25-S1 features a facility for CMN reduction; it has an external 
resistor in the power supply line common to all 128 channels, which collectively drives 
down their output in case of a CMN pulse. The remaining CMN can be determined and 
subsequently subtracted for each of these groups separately within the data of one event. 

During the intended operating time of ten years, the components of the CMS tracker 
will be exposed to an aggressive environment. The radiation level within the tracker 
volume will be extremely high, such that some components will encounter a fluence of up 
to 1.6 ■ 10^^ MeV-equivalent neutrons per cm^ [81] . As a consequence, the silicon has to 
be protected from radiation damage by cooling it to a temperature of -10°C. 

3.4 A Laser Test Facility for Silicon Microstrip Detectors 

To provide a stable and proper operation of the CMS tracker according to its performance 
specifications, care has to be taken that all assembled components fulfill strict quality 
requirements. Regarding the silicon microstrip detector modules, these requirements are 
summarized within a grading scheme as follows [H2]- Modules are classified as grade A 
if they have less than 1 % defective strips, as grade B if they have between 1 and 2 % 
defective strips, or grade C otherwise. Modules of grades A and B are used for tracker 
assembly with grade A preferred. In addition to the requirement regarding defective 
strips, modules may also be disqualified due to deviations in their electric properties, the 
response of the frontend electronics and the mechanical precision of assembly. 

3.4.1 Possible Module Defects 

After completing the ring 6 module assembly with the final steps of sensor to sensor and 
sensor to pitch adapter wirebonding, there are three main categories of flaws which may 
be responsible for the occurrence of a defective strip: 

• Open or missing wirebonds, which are caused by flawed calibration of the wire- 
bonding facilities or improper handling of the modules. In principle, this category 
also includes interruption of any signal lines, such as the pitch adapter conductors or 
the aluminum strips themselves. Open or missing wirebonds can largely be repaired. 

• Inter-strip short-circuits - ohmic contacts between adjacent strips - which al- 
most exclusively result from manufacturing defects. They degrade the detector res- 
olution and are in most cases irreparable defects. Hence, inter-strip short-circuits 
are considered in the module grading and must therefore be identified. 

• The so called pinholes, which represent ohmic contacts between the aluminum 
strips and the underlying p"*" implants through the insulating oxide/nitride layer. 
Pinholes may result from manufacturing defects or from flawed wirebonding. One 
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particular issue about pinholes is that already a single one allows current to flow 
from the sensor bulk into the frontend electronics, possibly leading to a breakdown 
of the corresponding 128 APV25-S1 channels at a time. Hence, strips affected by 
pinholes must be detected and detached from the hybrid by removing the wirebonds. 

A test facility for silicon microstrip detector modules must reliably detect defects from 
these three categories. It has to provide easy handling and, once the module is mounted 
for testing, fully automatic operation and high throughput are required. In total, 1061 
ring 6 modules (including spares) are produced for the SST, all of which are wirebonded 
at the Physikalisches Institut IB at the RWTH Aachen. A test facility for these modules, 
which meets the above requirements, has been developed and is described in the following. 

3.4.2 The Setup of the Module Test Facility 

The basic principle of the module test facility (MTF) is the generation of signals in the 
silicon sensors by laser illumination and their subsequent analysis. As will be shown, 
each defect leaves a signature in the signal of the affected strip, characteristic of each of 
the above categories and allowing their detection and tagging. Since the implants are 
capacitively coupled to the readout strips, the illuminating laser must be pulsed. 

The modules are read out via the ARcH system [771 ES]. ARC is a multicomponent 
readout system developed for module testing purposes and comprises the ARC controller 
board, two frontend hybrid adapters, an ISA card interface for personal computers, and 
the operating software package ARCS. The ARC board represents the core of the system 
and performs clock and trigger generation and distribution, data sampling and buffering, 
as well as slow control and voltage control for the IC on the hybrids. Two modules can 
be read out in parallel; their data are sampled by a group of three 8-bit analog-to-digital 
converters at optionally 20 or 40 MHz. The ARC board additionally accepts external 
trigger signals. Of the two frontend adapters, one houses amplifiers, line drivers and 
communication circuits while the other serves as a voltage regulator for the hybrids. The 
operating software ARCS is arranged in three levels: a hardware device driver, a package 
of communication routines written in the high-level programming language C++, and a 
graphical user interface. In the configuration implemented in the MTF, the ARC system 
can deliver data at a maximum rate of about 200 Hz. A detailed description of ARC is 
given in [77] . 

The electronic and optical signal lines of the MTF are depicted in Figure 13.91 They 
are grouped in two branches, a laser or signal generation branch and a branch for trigger 
generation and readout. Both branches start at the pulse generator's output jacks, which 
deliver two separate square-wave signals with a duty cycle of 50 % and a frequency of 
180 Hz. The first signal with positive voltage and an amplitude of 0.8 V is fed in the laser 
unit, which emits at a central wavelength of 1050.5 nm with an amplitude of 0.43 mW. 
This infrared light is guided via a single-mode optical fiber to the laser optics, which focus 
the beam on the silicon sensor surface at a distance of 50 mm. The resulting laser spot 
has a nearly Gaussian intensity profile with a minimum achievable width of 10 fim. For 

^ARC - APV readout controller 
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Fig. 3.9: Schematic view of the laser test facihty's optical and electrical circuitry. 



the purpose of module testing, the spot width has been fixed close to this minimum value, 
so that the strips can be targeted individually. 

Simultaneously, a second signal leaving the pulse generator with alternating voltage 
sign is fed into the trigger unit (gate generator), which transforms it into a series of stan- 
dardized logic NI1V0 signals with voltages of V and —0.8 V for logic and 1 respectively. 
The output is connected to the trigger input jack of the ARC board. On the arrival 
of a logic 1 pulse, the ARC board executes the readout sequence of the APV25-S1 at a 
constant phase with respect to the laser pulse inducing signals in the silicon sensor. The 
phase - i.e. the pipeline slot to be read out - is determined by the difference of the total 
delay times for signals traveling along the two branch lines and has been optimized by 
samphng tests. 

Fi gure 13 . 1 01 shows a photograph of the mechanical layout of the MTF. During the final 
steps of their assembly, the modules are mounted on standardized aluminum support 
plates for easy handling and interim storage. For laser testing, a module is placed in an 
aluminum/plastic composite box which is equipped with a slide- in slot for quick insertion 
of the support plate. The box itself is mounted on a rigid aluminum pole and carries the 
electrical connectors for the readout system and the depletion voltage supply. Across from 
the module housing box, a motorized linear stage is mounted horizontally, itself carrying 
a second vertically oriented motorized stage, to which the laser optics are attached. By 
these means, the optics can be moved in a plane parallel to the module and illuminate 
any point on the silicon sensors. Both stages are driven by servo motors and offer a 
positioning accuracy of less than 16 fim. The entire setup is housed in a lighttight cabinet 
and is mounted on an air-cushioned granite table to protect it from vibrations. Table 13.21 
summarizes the components of the MTF and their properties. 



3.4.3 The Test Procedure 

A personal computer (PC) is deployed for controlling the course of a module test session 
with the MTF. It addresses the two linear stages, operates and monitors the high voltage 
power supply, and stores the silicon module data from the ARC board. The custom 

''NIM ~ Nuclear Instrumentation Modules 
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Fig. 3.10: Photograph of the module test facihty. 



software package for the MTF is written in C++ and makes use of the ARCS hardware 
device driver and hybrid communication routine level, as well as a custom device driver 
for the linear stage controller board. Figure 13.111 illustrates the path followed by the laser 
spot in the plane of the module surface during the course of the fully automated testing 
procedure. In a first step, after insertion of the module into the MTF box, the depletion 
voltage for the module is slowly ramped up to 100 V to avoid high currents flowing into 
the capacitance of the growing depletion zones around the p-n junctions. The ramp time 
is set to 30 s. Subsequently, the laser spot is moved to a position away from the silicon 
sensor surfaces in order to take pedestal data without laser illumination. 

The pedestal Pi of a sensitive detector element i is defined as its average output level if 
no signal is present. Due to manufacturing tolerances of the hybrid IC, the pedestals of 
the individual strips differ slightly from each other and must therefore be determined for 
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Fig. 3.11: Schematic view of the path which the laser spot follows during the test procedure 
in the plane of the module surface (solid blue line). 



later subtraction from the signal data. For this purpose, the module is read out a certain 
amount of times n with the laser in the offside position and the resulting data Di^n are 
stored for each strip i. The pedestals are then calculated according to 

1 " 

= -5]A,fc. (3.6) 

fc=i 



controller 


Intel P3 866 MHz personal computer 


linear stages 


OWIS Limes 90 (vertical) and Limes 120 (horizontal) DC servo stages 
with 150 mm and 200 mm travel, respectively, and <16//m position- 
ing accuracy; OWIS DSC 1000 servo amplifier actuated by a Galil 
DMC 1040 motion controller for ISA bus 


pulse generator 


TTi TG550 function generator 


laser source 


Custom built single- mode diode laser unit, center wavelength A = 
1050.5 nm, with output power < ImW and Ins rise time 


trigger unit 


LeCroy model 222 gate generator 


laser optics 


Schafter+Kirchhoff micro focus (/ = 50 mm) mounted on a fiber 
collimator (/ = 18.4mm); 1.5m long single-mode fiber with 7.5 /im 
core diameter; 8° angled polish connectors 


module readout 


ARC system including ARC board, frontend adapters, ISA computer 
interface and cabling 


depletion voltage 
supply 


ISEG SHQ 222m power supply with RS 232 interface 


software 


Linux kernel 2.2.14 operating system, custom driver for linear stages, 
ARCS driver, custom software for module test operation 



Tab. 3.2: Details of the components of the laser test facility. 
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In the case of the MTF, a value of n = 10000 has been chosen to determine the pedestals 
with high accuracy. 

Since at the beginning of a test session the state and position of the linear stages must 
be assumed to be unknown, an absolute position calibration is performed as the next step. 
For this, the two stages are moved towards the upper (vertical stage) or right (horizontal 
stage) boundary of their range, until they are stopped by their internal limit switches. 
This position is reproducible within the 16 /xm accuracy of the servomotors and is then 
defined as the reference point by resetting the position encoders to zero. 

After calibration, the MTF is ready for laser testing. The laser spot is moved succes- 
sively over both sensors in a direction perpendicular to the strips at a constant speed of 
200 /im/s. During the scans the individual strip signals are read out and stored with the 
APV25-S1 operated in peak mode. To allow checking the full strip length for possible 
defects, the sensors are scanned at the end far from the readout electronics. Scanning 
both sensors separately gives additional information about the identity of the sensor or 
wirebond which is affected by a particular defect. The channel numbering used in the 
following is indicated on the left side of Figure 13.111 and starts at zero. The laser moves 
along the direction of increasing strip numbers in both cases. The final output of the 
procedure is an array of 512 values for each scan - the set of the maximum signal deliv- 
ered by each of the 512 strips. From these numbers, the corresponding pedestal values 
are subtracted. Including pedestal data collection and calibration of the linear stages, the 
total time needed for a test session with the MTF amounts to approximately 900 seconds 
per ring 6 module. 



3.4.4 Investigation of Artificially Bonded Defects 

To investigate the performance of the MTF and to assure its applicability for reliable 
detection of module defects, a reference module has been prepared with artificially bonded 
defects. Two pairs of short-circuited strips were simulated by placing a wirebond between 
the AC pads of the strip pair. Additionally, three pinholes resulted by wirebonding the 
AC and DC pads of a strip, thus establishing a perfect short-circuit between the surface 
aluminum and the underlying implant. The wirebonding was performed at the far end 
of sensor W6A, at maximum distance from the hybrid. Furthermore, one open wirebond 
between the sensor W6B and the pitch adapter appears on the reference module. Table EH] 
summarizes the artificial defects and gives the number of the affected strips on the module. 



defect type 


number 


strip numbers 


inter-strip short-circuits 


2 


269-270, 288-289 


pinholes 


3 


309, 329, 349 


open wirebond 


1 


411 



Tab. 3.3 : The type and number of artificially bonded defects on the reference module and the 
affected strip numbers. 
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Prior to the laser testing procedure, the pedestals and the CMN corrected noise of 
the reference module have been determined with the ARCS software [7S]. Since the 
analog inputs of the APV25-S1 are bundled in four blocks of 32 channels each, the CMN 
correction algorithms implemented in ARCS are based on this group definition [77] . The 
results are depicted in Figure 13.121 The left part of the figure shows the strip pedestals, 
which exhibit a typical structure with sharp edges at the APV25-S1 borders and a slope 
within the range of the individual APV25-S1, which is probably due to a drop of supply 
voltage across the channels. No hint for any defective strip can be observed in these data. 




Fig. 3.12: The pedestals (left) and the common mode subtracted noise for the two APV25- 
Sl readout modes (light) for the MTF reference module. The APV25-S1 border channels are 
indicated by red dashed marks. 



The right side of Figure 13.121 displays the noise of the strips of the reference module 
for the two APV25-S1 readout modes. In both cases, the noise data show a largely flat 
common baseline with a series of spikes of varying amplitude in both directions and a lower 
plateau in the strip range of one APV25-S1. The artificially bonded defects, indicated by 
the corresponding strip numbers, are clearly visible as distinct peaks or dips and can thus 
in principle be tagged as faults. However, a number of effects complicate the identification 
of defective strips in the noise data. 

First, the channels on the APV25-S1 borders and especially the outermost strips close to 
the edge of the module show an increased noise level. This behavior is mostly due to their 
position which affects the capacitive couplings to their neighbors and other circuit paths 
on the module (e.g. to the guard or bias ring, see Figure [H^ . Despite the peculiarities 
shown by border and edge strips in the noise data, they are not to be considered as 
defective, thus making it necessary to distinguish them from the truly faulty channels. 
Furthermore, the signature of the different categories of defects is not unique in the noise 
data. Obviously, as apparent for example in the behavior of the pair of short-circuited 
strips 288-289 and the open wirebond (number 411), the signature varies with the readout 
mode, as the channels exhibit a lower or higher than average noise level in the two cases, 
respectively. The reason for this behavior lies in the CMN correction facility of the 
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APV25-S1, which apphes also to channels with lower noise as for example those affected 
by an open wirebond. The correction may then lead to a higher apparent noise level, 
depending on the readout mode. In addition, the presence of multiple pinholes in the 
strip range of a single APV25-S1 leads to a significant lowering of the noise baseline due 
to an increased current flow between the silicon and the APV25-S1, which complicates 
the detection of defects. Finally, the signatures of defects in the noise data are highly 
sensitive to the environmental circumstances under which the module is operated, such 
as grounding and electromagnetic shielding, making them generally unpredictable and 
unreliable for testing purposes. 
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Fig. 3.13: Laser scan results from sensor W6A of the reference module (left); enlarged view 
of the strip interval exhibiting conspicuous strip signals from artificially bonded defects. The 
strip numbers are indicated (right). 



Figure 13.131 displays the result of the laser testing procedure as it was applied to sensor 
W6A of the reference module. The overview histogram on the left side contains the 
maximum laser scan signal for each of the 512 channels; the right side shows an enlarged 
view of the strip range with bonded defects. As apparent in the overview, the normal 
channels show a nearly constant response to the laser illumination of about 110 adc counts. 
However, a slightly negative slope remains, which is caused by interference effects from 
light reflected at the metalized backplane of the sensors. Such effects are due to slight 
variations of the sensor thickness of the order of a few micrometers and have already 
been observed previously [HI]. AH artificially bonded defective channels clearly appear as 
distinct negative spikes from the common baseline and consistently show a lower response, 
independent of the readout mode. An indication for the type of defect can be obtained 
from the laser scan signature in the following way. 

• Short-circuited strips are recognizable through their pairwise occurrence. Besides, 
since they are connected, the charge induced on one strip is shared among the two 
readout channels and leads to a reduced signal response for both. As a consequence 
of the increased strip capacitance at each channel input, the noise level is higher. 
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so that the resulting laser response of two particular pairs of short-circuited strips 
may significantly differ. 

• Pinholes represent direct connections from the aluminum strips to the underlying 
p"*" implants, permanently forcing the channel inputs of the affected strips to the 
implant potential. Therefore, no signal is observed at any time and the noise level is 
very small. Pinholes are thus identifiable through their particularly small response 
to laser illumination. 

• Open wirebonds, similar to pinholes, prohibit the observation of any signal in re- 
sponse to laser illumination. However, due to the noisy behavior of the affected 
channels, resulting from the CMN correction described above, their mean observed 
response is typically higher than in the case of a pinhole. 

It must be pointed out that, as a consequence of the statistical character of signal 
noise, the above criteria for identification of the type of defect must not be considered to 
be unambiguous in every case. The main purpose of the MTF is reliable defect detection. 
Figure 13.141 displays the distribution of the response for the channels of the reference 
module. The bulk of flawless strips show an approximately Gaussian distribution with a 
mean of 107 and a standard deviation of 4.3 adc counts. A small tail on the right side is 
caused by the slightly higher response of the module's edge strips. The region populated 
by the defective channels on the left side is clearly separated from the bulk by more than 8 
standard deviations. Hence, all artificially bonded defects are reliably identifiable through 
their low response to laser illumination. 



3.4.5 Results from the TEC W6 Prototype Series 

The TEC collaboration produced a prototype series - the express-line 2 - of nine ring 6 
modules. The main purpose of the express-line modules was to allow the development, im- 
provement and validation of module test procedures prior to series production. Contrary 
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Fig. 3.14 '- Distribution of the laser response of the reference module channels. The red dashed 
line represents a Gaussian fit to the bulk of faultless strips, whose parameters are given in the 
legend. Defective channels are indicated by their strip numbers. 
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to the final module design, the express-line modules are equipped with ceramic hybrids 
instead of the kapton/copper ones, which, however, has no effect on their readout prop- 
erties in this procedure. Individual modules are distinguishable through identification 
numbers on the support frame, which are listed in Table I3.4[ 

All nine modules of the express-line have been tested with the MTF according to the 
procedure described in § 13.4. 31 The results of the individual scans of the W6A sensors are 
depicted in Figure 13.151 Confirming the earlier results on the reference module, defective 
strips are clearly visible through their strongly reduced response to laser light. Four of 
the modules turned out to have no faulty channels at all; their laser scan results are 
summarized in panel a) of the figure. Three of the remaining modules (M2, M3 and M6) 
have one defective channel each, while modules M4 and M7 each have two faulty channels. 
According to the grading scheme introduced in § 13.41 based on the number of defective 
strips, all modules must be considered as grade A. Nevertheless, modules M5, M6 and 
M7 were later classified as grade C due to their high leakage current [53] . 



module 


frame id number 


faulty channels 


Ml 


30200020000501 




M2 


30200020000502 




M3 


30200020000641 


242^ 


M4 


30200020000643 


447P, 509 P 


M5 


30200020000503 




M6 


30200020000644 


357P 


M7 


30200020000640 


77P, 89^ 


M8 


30200020000642 




M9 


30200020000645 





Tab. 3.4 '■ The results of the express-Hne module tests with the MTF. Faulty channels are 
indicated by the type of defect: p for a pinhole and ^ for an open wirebond. 



Table 13.41 gives an overview of the modules and the defective channels detected. From 
the seven defects listed, two have been identified as open wirebonds by visual inspection, 
and four represent pinholes. The failure of channel 242 on M3 is probably due to a 
broken pitch adapter line [85]. No short-circuited strips were observed. The results of the 
express-line module tests with the MTF turn out to be consistent with those from other 
facilities using different testing techniques [77] . 

Following the expess-line module tests the MTF is established as a practicable and 
reliable facility for the detection of defective strips on the TEC modules. Furthermore, 
with its high positioning accuracy and stability of the laser intensity, the MTF is well 
suited as a precision tool for the in-depth investigation of silicon microstrip sensor prop- 
erties. However, the disadvantage of the facility are the high costs of the hardware. In 
order to provide all laboratories involved in module testing with identical test facilities, 
the decision was made to use an extension to the ARC system [7S] for testing during the 
series production of the modules. Besides test procedures based on measurements of the 
sensor leakage current, APV25-S1 amplifier gain linearity, noise distributions and fron- 
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tend electronics performance, the extended ARC can perform illumination of the silicon 
sensors with light-emitting diodes. 

The assembly of the CMS Tracker Endcaps was finished in November 2006 and first 
operation of the LHC is scheduled for 2007. The following chapters describe the search 
for signatures of supersymmetry in the cosmic ray spectra and the analysis of data from 
the AMS-01 experiment, which bases on the reconstruction of multi-particle events in a 
large-scale application of a silicon microstrip detector which was operated in space. The 
subject is introduced with a discussion of cosmic rays in the next chapter. 
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Fig. 3.15: Laser scan results from W6A sensors of the express-line modules. Display a) 
combines the scans of the four modules without defects Ml, M5, M8 and M9. For the sake of 
clarity, the common baselines of the channels have been shifted downward by 120, 80, 40 or 
adc counts, respectively. The remaining histograms b) to f) show the individual results for 
modules M2, M3, M4, M6 and M7, with the defective strips indicated by their strip numbers. 
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The Earth is continuously exposed to a flux of particles of extraterrestric origin, referred 
to as cosmic rays. This was discovered in 1912 by V. Hess, who observed an increasing dis- 
charge rate of electrometers towards higher altitude during a series of balloon flights [HS] . 
Historically, the investigation of cosmic rays has played an important role in the devel- 
opment of particle physics, as they were the only source of high energy particles before 
the emergence of accelerators. Thus, the positron [H7j, the muon [HH], the pion [HH], the 
kaon and other strange mesons have been discovered in cosmic rays. 

Experimental knowledge about cosmic rays spans a very wide energy range. Particles 
from the solar wind typically propagate with a few keV, while, on the other hand, cosmic 
rays have been observed with energies higher than 10^'^ eV [9l] . Moreover, the abundance 
of cosmic rays strongly varies as a function of energy. Over large energy intervals the 
spectra follow a steep power law with a spectral index of roughly 3, meaning that the 
intensity of cosmic rays decreases by a factor of 1000 for each decade in energy |95j . 

Regarding the search for dark matter, one of the most important properties of the 
spectral shapes of cosmic rays is their smoothness on small scales. As will be shown in 
this chapter, the cosmic ray background produced in conventional high energy processes 
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Fig. - Total cosmic ray flux as a function of energy [9T] (left); fluxes of some of the most 
abundant cosmic ray components as measured by AMS-01 WI\ and HEAT-e^ [3] (right). 
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in the universe is expected to lack any structure. Reaction chains involving dark matter 
may be additional sources of standard model particles and would cause an excess within 
limited energy ranges. Of particular interest are the fluxes cosmic ray particles of low 
abundance, such as the positron, since exotic sources would lead to significant excesses. 



4.1 Cosmic Ray Composition and Spectra 

The particles reaching the top of the Earth's atmosphere are mainly electrons and fully 
ionized light atomic nuclei, such as hydrogen and helium. Nuclei with higher charge and 
other particles amount to less than 1 % of the total intensity. Table 14.11 gives an overview 
of the approximate abundance of the most dominant cosmic ray components at 1 GeV 
energy. Antinuclei have not been observed yet The chemical composition of cosmic 
rays is similar to that of the sun, which indicates that they are predominantly of stellar 
origin. Exceptions are the nuclei of elements from Li to B and from Sc to Mn, which are 
overabundant as the products of spallation of carbon and oxygen and of iron, respectively. 



Component 


Abundance 


[%] 


Component 


Abundance 


[%] 


Hydrogen 


90 


Positrons 


0.4 


Helium 


4.6 


Carbon 


0.2 


Electrons 


2.3 


Oxygen 


0.2 


Deuterons 


1.7 


Other nuclei 


0.3 



Tab. - Approximate abundances of the most dominant components of cosmic rays at 
1 GeV (or IGeV/nucleon for nuclei) energy, in percent of the total intensity. Based on data 

from [adellMllMlETlEH]. 



Due to the deflection of charged particles in the galactic and solar magnetic fields, the 
arrival directions of cosmic rays are randomized at detection, so that the fluxes are widely 
isotropic. Their energy spectra above 10 GeV can be described by a segmented power law 
of the form 

dN 

with the spectral index 7 as an essential parameter. Up to energies of 10^^ eV 7 is equal 
to 2.7 and then steepens to 3.0 for higher energies (see Figure SH], left). The break in the 
spectra at 10^^ eV is commonly referred to as the knee. At approximately 10^^ eV - the 
ankle - and above, the spectra seem to flatten again. This behavior is independent of the 
incident direction, i.e. the spectral shapes are also isotropic. Below 10 GeV, the spectra 
are modulated by the solar and geomagnetic influence, which is the subject of § 14.31 

The reason for the two breaks in the spectra has not been fully understood as yet [95] . 
The conventional explanation is that they represent the power limits of different accelera- 
tion mechanisms at successive energy scales. Obviously, these mechanisms must be based 
on a common principle, since the spectral shape is universal over large energy ranges. 
It was proposed in 1949 by E. Fermi [99], that the spectral shape may originate from 
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repeated finite energy gains of particles in collisions with moving objects, e.g. magnetic 
fields in shock waves of supernovae. Let a particle of initial energy Eq undergo such a 
repeated acceleration and gain a fraction $,E oi its energy E each time. Then, after n 
collisions, it will have 

E„ = Eo(l+0", (4.2) 

giving for n 

n = \n{E/Eo)/Hl + 0- (4.3) 

Assuming a probability for a particle to escape the acceleration mechanism, the number 
of particles with an energy larger than E is 

Ni> i?) oc E (1 - = ^ • (4.4) 

m=n Pw 

Substituting n from eq. 14.31 yields 

N(>E)^ — (^) \a^p^/i, (4.5) 

which leads to the observed form of the spectra |100j . This is the principle on which the 
so called second order Fermi acceleration mechanism is based. However, in this simple 
form, it cannot sufficiently explain the isotropy of the observed spectra. Furthermore, 
it requires large initial energies and is very inefficient, so that cosmic rays of highest 
energy are unlikely to have been accelerated in this way. Consequently, Fermi's theory 
has undergone fundamental further development in the past decades |1U1] . 

Figure HIT] (right) shows the spectra of protons, helium nuclei, electrons and positrons 
with energies up to approximately 100 GeV. It is apparent that in this energy regime 
protons are by far the most abundant cosmic ray component. In particular, their flux 
exceeds that of positrons by a factor of 10^ — 10^. The ratio of the positron to electron 
flux varies from roughly 1:7 to 1:20. Moreover, the energy spectra of the leptons (7 ~ 3.4) 
are slightly steeper than those of nuclei (7 ^ 2.7). 



4.2 Origin and Propagation of Cosmic Rays 

Cosmic rays can be classified in two main categories of different origin, solar and galactic 
cosmic rays. As obvious from these terms, solar cosmic rays have their origin and are 
accelerated in the vicinity of the sun, while galactic cosmic rays are produced somewhere 
in the galactic volume. However, it is often assumed that cosmic rays of highest energies 
may even originate from extragalactic sources. On their way to Earth, cosmic ray particles 
are involved in a plenitude of processes, including elastic and inelastic scattering off the 
interstellar medium, deflection by magnetic flelds and deceleration through emission of 
radiation. 

Solar cosmic rays mostly consist of ionized coronal material which is accelerated by 
shock waves associated with the solar wind |102] . The abundance of these particles de- 
creases steeply with energy. At low energies, below 100 MeV, they constitute the major 
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part of the observable cosmic ray flux, but in the GeV energy range their persistent 
contribution is negligible. Only in the rare case of intense solar flares may particles be 
accelerated up to 1 GeV or even more [1U3] . These events have typical durations of the 
order of minutes or hours. The composition of cosmic rays produced in solar flares is 
quite similar to that of galactic cosmic rays [104] . 

Supernovae provide the dominant power input to the galaxy. An average supernova 
explosion releases between 1 and 10 solar masses of material and a total energy of 10^^ J. 
Taking into account a supernova explosion frequency in the galaxy of roughly 3x 10^ Ma~^, 
it follows that each supernova would only have to contribute ~ 3 % of its output power to 
the total cosmic ray powei0 of 10^^ W [104] . For this reason supernovae are believed to be 
the dominant accelerators of cosmic rays to energies up to 10^^ eV. Since particles are scat- 
tered by interstellar gas and therefore their energy declines inversely with the cube root of 
the number density, they cannot be accelerated too soon after the explosion. Rather than 
that, the acceleration takes place when the interstellar gas is passed by the shock waves 
emerging from the supernova. Evidence comes from obervations of supernova remnants 
in the keV X-ray band |105] . yielding a flat spectrum which can be interpreted in terms 
of bremsstrahlung emission from protons or electrons. Furthermore, recent results from 
imaging atmospheric Cerenkov telescopes reveal that these remnants also emit gamma 
rays with up to TeV energies |106] , in agreement with the X-ray measurements regarding 
the spatial location of the sources. Besides supernovae, other candidates for sources of 
cosmic rays are pulsars, close binary systems and stellar winds |107j . Particles which are 
produced and accelerated in the vicinity of these sources are, with respect to their origin, 
commonly referred to as primary cosmic rays. 

Between its acceleration and its detection a galactic cosmic ray particle traverses an ave- 
rage amount of interstellar matter of lOgcm"^ |108] within typically 20 Ma of time [109] . 
The propagation volume is permeated by the galactic magnetic field, whose regular (large 
scale) component has a strength in the order of 5 /iG |110] , however with turbulent local 
variations. As a consequence, in addition to their convective thermodynamic motion, 
charged particles diffuse due to interactions with the magnetic field inhomogeneities. 
While propagating through the galaxy, cosmic rays also lose energy. For nuclei, the 
dominant energy loss processes are ionization in the interstellar medium and scatter- 
ing off thermal electrons in plasmas |109] . Electrons and positrons additionally undergo 
bremsstrahlung, synchrotron radiation and inverse Compton scattering with starlight or 
with the cosmic microwave background [111] . Assuming the presence of a galactic wind, 
i.e. a constant flow of magnetic irregularities, particles may also gain energy while pro- 
pagating from shock acceleration (diffuse reacceleration) . 

The destruction of primary nuclei via spallation in the interstellar medium gives rise 
to secondary nuclei and rare isotopes. Photons and antiparticles are also rare in cosmic 
rays. Apparently, they are not directly produced in the cosmic ray sources, but emerge 
as secondary products from processes involving primary cosmic rays, and are thus called 
secondary cosmic rays (see Figure 1121 left). Antiprotons largely stem from proton-proton, 
proton-nucleus and nucleus- nucleus interactions |112j . Photons are decay products of 
neutral pions from such processes and also result from synchrotron radiation of high 



i.e. the estimated total power leaving the galactic disk in form of cosmic rays 
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Fig. 4-2: Principle of primary and secondary cosmic ray production (left); model calculation 
of the positron fraction /{e'^ + for purely secondary positron production from [lllj (solid 
line) together with recent data from AMS-01 g], HEAT-e^ and HEAT-pbar [7] (right). 



energetic particles or inverse Compton scattering. Positrons are mainly created by protons 
in reactions of the type pp ^ + X or pp ^ + X, where the pions and kaons decay 
into positrons via muons. Furthermore, the 7e"'"e~ decay channel of neutral pions from 
reactions of the type pp ^ + X contributes to the positron flux. In addition to protons, 
light nuclei may also be involved in these interactions [113] . 

The cross sections of the processes which contribute to the positron flux are strongly 
energy dependent and the spectra of their individual contributions are smeared by energy 
losses during particle propagation. The ratio of the positron and electron fluxes, commonly 
written in terms of the positron fraction e^/(e^ + e~), is thus a smooth function of 
momentum. Figure l4]2] (right ) shows a model calculation of the positron fraction for purely 
secondary positron production (without reacceleration) together with recent experimental 
data. In the range from a few hundred MeV/c to 50GeV/c it decreases monotonously 
from 0.2 to roughly 0.03 and then flattens out. In particular, for purely secondary positron 
production, the positron fraction does not exhibit any small scale features, peaks or dips. 



4.3 Cosmic Rays in the Earth's Vicinity 

Within the heliosphere, cosmic rays are influenced by the solar wind, a continuous stream 
of charged particles emanating from the sun. The magnetic field associated with this 
stream has irregularities which scatter particles in a frame of reference moving outwards 
of the solar system. As a result, cosmic rays are effectively decelerated (or even deflected) 
in the vicinity of the sun, an effect which is referred to as solar modulation. According 
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to the so called force field approximation [114] . the influence of solar modulation can be 
parameterized with a single parameter 0, which represents an effective potential and is 
given in units of volts. Thus a particle with charge Ze experiences an energy loss Ze(j) 
when approaching the Earth from outside the heliosphere. As displayed in Figure 14.31 
the solar modulation parameter is time dependent and follows the 11 year solar cycle, 
typically varying from 300 MV to 1300 MV. As can be seen e.g. in Figure W7L\ (right), the 
solar modulation is apparent as a decrease with respect to the power law expectation of the 
cosmic ray fluxes below approximately 10 GeV. There is evidence that solar modulation 
is charge sign dependent at low particle energies [115] . 
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Fig. 4-3 '■ The solar modulation parameter cj) (solid line) with 68% confidence intervals (grey 
area) as a function of time |116j . The circles denote reference data. 



Immersed in the heliospheric particle streams is the magnetosphere of the Earth. The 
rather complex magnetic field generated by the Earth's body (internal field) is commonly 
treated as the derivative of a scalar potential, B = — grad V, with V expanded in terms 
of spherical harmonics, 

V = aJ2Y.i9n cos mcp + sin m0) - (cos 9) . (4.6) 

n=l m.=0 

Here, a is the mean radius of the Earth, {r,(j),6) are geocentric coordinates, P^{cos6) 
are the associated Legendre functions, and (5'JJ\ h^) are the Gauss coefficients describing 
the field contribution of the separate terms in the sum [117] . Additional terms may 
also incorporate external magnetic field sources, such as particle streams generated by 
interactions with the solar wind. The dominant terms in eq. (14. 6 p are related to n = 1, 
i.e. the geomagnetic field is in first order that of a dipole, tilted by about 168.5° with 
respect to the rotation axis of the EartlS and displaced approximately 400 km from its 
center. The field strength on the Earth's surface varies between 20 fiT and 70 fiT. As 
shown in Figure 14.41 (left), at distances of a few Earth radii, the geomagnetic field is 

^Thc geomagnetic north pole lies in the southern geographic hemisphere. 
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distorted by the solar wind, a highly dynamic process which gives rise to a shock front 
(bow shock) along the border between interplanetary magnetic fields and the Earth's 
magnetosphere (magnetopause) . 




Fig. 4-4 ■ Schematic view of the Earth's magnetosphere and the van Ahen belts, as distorted 
by the solar wind (left); trajectory of a charged particle in the geomagnetic field (right) |118j . 

Low energy charged particles can be trapped in the geomagnetic field. The path which 
such a particle travels is schematically illustrated in Figure (right). The basic tra- 
jectory is a helix following the field lines, with an increasing radius of curvature of its 
gyromotion {gyroradius) as the particle moves further away from the Earth. Near the 
Earth's surface, where the field strength is larger, the gyroradius decreases. The parti- 
cle's energy then becomes dominated by the gyromotion contribution at the expense of 
the forward motion along the field line. Consequently, the particle is defiected back along 
its path, only to repeat the process at the other end. As the trapped particle bounces 
back and forth along the magnetic field, it also drifts around the Earth because, in part, 
its gyroradius is larger on the outer part of its trajectory than on the inner part. This 
effect causes negatively charged particles to drift in the direction of the Earth's rotation 
(eastward) and positively charged ones to drift in the opposite direction |118j . 

It follows that, in the dipole approximation, the region which can be occupied by 
trapped particles forms a nearly toroid-shaped volume symmetrical about the dipole axis. 
The Van Allen radiation belts surround the Earth at a distance of roughly 1.5 and 4-5 
Earth radii from its center [119j . While the inner belt mostly consists of protons with 
energies larger than 100 MeV and is relatively stable, the outer belt is largely populated 
with electrons ranging in energy from 400 keV to above 15 MeV, and dynamically de- 
formed by the solar wind. Protons from the inner belt can interact with atmospheric 
molecules, giving rise to secondary particles (atmospheric secondaries) such as, amongst 
others, positrons. 

The ability of a cosmic ray particle to penetrate the Earth's magnetosphere is uniquely 
determined by its rigidity, the momentum divided by the charge, given in units of volts. 
Cosmic ray particles with rigidities below a certain geomagnetic cutoff are deflected by 
the geomagnetic fleld and cannot reach the vicinity of the Earth. The value of this cutoff 
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Longitude [de^] 

Fig. 4-5 : The vertical cutoff rigidity at an altitude of 450 km under magnetically quiet condi- 
tions. Based on data from |120j . 



not only depends on the geomagnetic coordinates of the particle's entry point into the 
magnetosphere, but also on its incident direction. 

Figure 113] shows the cutoff rigidity for vertical incidence at an altitude of 450 km above 
the Earth's surface under magnetically quiet conditions |120] . It ranges from a few MV 
in the polar regions up to approximately 15 GV near the equator. Similar to the solar 
modulation effect, the geomagnetic cutoff modulates the spectra of cosmic rays having 
rigidities in this range. Particles detected in an experiment with rigidities below their 
cutoff cannot be of galactic origin, but must stem from within the magnetosphere. They 
are mostly produced as atmospheric secondaries and are trapped in the Earth's radiation 
belts. 

4.4 Signatures of SUSY Dark Matter in the Cosmic Ray 
Spectra 

As stated already in sec. 12.41 neutralinos as Majorana particles could annihilate pairwise 
into SM particles, which then propagate as cosmic rays and are thus accessible to ex- 
periments. This process would constitute an additional primary source of cosmic rays 
with a unique spectral shape, and may thus appear as an excess of particles over con- 
ventional expectations. However, since the corresponding source strength is assumed to 
be small compared to that of supernovae, the dark matter annihilation signal should be 
sufficiently significant only in the spectra of cosmic ray particles which are exclusively pro- 
duced through secondary processes at low rates, such as positrons, antiprotons or gamma 
rays (see § l4.2p . 

Figure 14.61 shows the Feynman diagrams for the dominant neutralino annihilation chan- 
nels. The process can proceed through the t-channel exchange of sfermions, charginos or 
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Fig. 4-6 '■ The dominant annihilation channels for the lightest neutralino x |47j . 



neutralinos, or through s-channel annihilation via Z° or pseudoscalar Higgs bosons. The 
final states are either fermions or gauge bosons. The total cross section for neutralino an- 
nihilation strongly depends on the individual channel amplitudes, and thus on the choice 
of supersymmetric parameters. For example, in the case of fermion final state diagrams, 
the amplitude is proportional to the fermion mass in the low velocity limit and inversely 
proportional to the squared mass of the sfermion /. Furthermore, the amplitude for Higgs 
exchange is proportional to tan j3 for down-type quarks [17] . This implies that at values 
of tan/5 > 5 neutralino annihilation is dominated by s-channel Higgs exchange with bb 
pairs as the final state. 

Subsequent to their generation, 6-quarks hadronize and give rise to decay cascades 
in which, amongst others, positrons and antiprotons are produced, as well as typically 
30-40 photons per annihilation [122] . Additional positrons originate from the decay of 
gauge bosons and of r and fi. Direct annihilation into electron positron pairs is helicity 
suppressed [13|. As an example. Figure |4?7I (leR) shows the raw spectra of positrons (prior 
to their propagation through the interstellar medium) from the different annihilation 
modes of neutralinos with a mass of 300 GeV/c^. The spectral shapes differ significantly. 




Fig. 4-'^- The positron spectrum from neutralino annihilations, prior to propagation, for 
selected annihilation modes and a neutralino mass of 300GeV/c2 [M] (left); positron fraction 
data and simulation results for purely secondary positron production (background) and with 
additional positrons from neutralino annihilation for two mSUGRA parameter sets, resulting in 
either a higgsino- or wino-like neutralino [121j (right). 
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since annihilations to heavy quarks and gauge bosons produce a much softer spectrum 
than the r or /i channels. The shape of the overall resulting spectrum is a mixture of the 
above and is determined by the cross section of the individual modes, which depend on the 
supersymmetric parameters. The positrons then propagate through the galactic volume 
and experience energy loss and deflection due to scattering and interactions with magnetic 
fields, which finally result in a smearing of the initial raw spectrum when observed near 
Earth. 

Quark hadronization functions and branching ratios of the decay modes involved in 
positron production from neutralino annihilations are known to a high degree of preci- 
sion [32] • Using models of particle propagation in the interstellar medium, the spectra of 
secondary particles near Earth can be simulated and compared to flux measurements. Fig- 
ure HTl (I'ight) displays the positron fraction as a function of the particle energy together 
with data from HEAT-e^ [3] and HEAT-pbar |123] . When compared to the background 
expectation, the data show an indication for an excess number of positrons at energies 
above 6 GeV. The solid and dashed lines represent simulations of the positron fraction 
including contributions from neutralino annihilation for two particular sets of supersym- 
metric parameters |121j . The spectral shape differs mainly due to the choice of neutralino 
mass. In the one case, the neutralino with a mass of 91 GeV/c^ is Higgsino-like, meaning 
that in the linear combination of eq. fl2.7p the coefficients C31 and C41 are large compared 
to the bino and wino contributions. In the other scenario, the 131GeV/c^ neutralino is 
wino-like with a dominant contribution from C21. In contrast to the bino-like neutralino, 
which largely annihilates to heavy fermions, gauge boson final states are preferred in the 
two cases above [12]. 

Both results involving dark matter annihilation clearly improve the agreement with the 
data compared to the background-only prediction. Clearly, in order to further constrain 




Fig. Data of the diffuse gamma ray energy spectrum fitted with the simulated background 
from different processes and an annihilation signal from a 50-70 GeV /c^ neutralino (left) j47j ; 
same for antiprotons with an assumed neutralino mass of 207GeV/c^ (right) |124j . 
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the mass range of the neutrahno, better data is needed, especially at higher energies in 
order to capture the trailing edge of the dark matter signal. The simulated positron flux 
resulting from annihilation is almost always too small to produce a visible signal. However, 
it may be increased by a boost factor B, which can be explained by assuming a clumpy 
nature of the dark matter distribution. In particular, since the neutralino annihilation 
rate is proportional to the square of its density, a signal enhancement is expected in high 
density regions like the center of the Galaxy. The boost factor is thus proportional to the 
mean deviation of the squared density p, i? oc (p^) / (p)^. 

Figure 14.81 shows the simulated background spectra of gamma rays and antiprotons 
compared to experimental results. In both cases the background alone turns out to be 
inconsistent with the measured fluxes, so that an additional source of particles is assumed 
to fit the data. In the gamma ray spectrum on the left side of Figure ITSl a signal from the 
annihilation of neutralinos with masses of 50-70 GeV/c^ is found to explain the observed 
excess well [17]. The right side of Figure US] shows the measured spectrum of antiprotons, 
which can again be fitted with the assumption of neutralino annihilation, however with a 
significantly different neutralino mass of 207GeV/c^ |124] . 

Despite the hints for the presence of a dark matter signal in the spectra of secondary 
particles, the neutralino hypothesis remains controversial. It has been pointed out re- 
cently that the excess of gamma rays, if explained as a signal of neutralino annihilation 
as shown in Figure I4l8] (left), may not be compatible with the observed abundance of 
antiprotons |125j . Other possible contributions to the cosmic ray positron flux have been 
proposed, such as electron-positron pairs from annihilations of hypothetical Kaluza-Klein 
particle states [12], or from conversion of synchroton photons emitted by galactic pul- 
sars |126] . Consequently, more data is necessary to allow further development of these 
models. While new experiments are yet under construction, one possible approach is 
the reanalysis of data from past experiments such as AMS-01 with new techniques. The 
description of the AMS-01 detector is subject of the next chapter. 
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5.1 AMS-01 Detector Construction 

As a prototype for the AMS-02 experiment, the AMS-01 detector was flown on board 
the space shuttle Discovery in a near Earth orbit during the STS-91 mission from June 
2nd to 12th, 1998. Its main purpose was to verify the feasibihty of engineering a space 
borne high-energy particle detector, launching it into orbit, and operating it safely under 
space environment conditions. Furthermore, the main design principles of the future 
AMS-02 experiment were tested. During its flight, AMS-01 has recorded a large amount 
of data for the determination of particle fluxes in the Earth's vicinity, leading to precise 
proton [HI [92] and helium [93] flux spectra, measurements of the characteristics of leptons 
in near earth orbit [2], and to a significant improvement of the upper limit on the relative 
flux of antihelium to helium in cosmic rays. 




Fig. 5.1: The AMS-01 detector: schematic cross section in the y-z plane (left) and three 
dimensional illustration (right). 

The structure of AMS-01 is illustrated in Figure 15.11 Six layers of double sided sil- 
icon strip sensors enclosed in a cylindrical permanent magnet formed the experiment's 
tracking chamber. The particle flight direction and velocity was measured with plastic 
scintillator paddles arranged in two double planes, one above and one below the mag- 
net (time of flight hodoscopes, TOF). Additionally, the inner surface of the magnet was 
covered with scintillators to reject background from particles which traverse the magnet 
wall or the support structure (anticoincidence or veto counters). As a complement to 
the velocity measurement with the hodoscopes, two layers of threshold Cerenkov coun- 
ters were situated at the lower end of the detector. In order to keep very low energy 
particles (up to several MeV) from entering, AMS-01 was shielded by thin carbon fiber 
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walls with a thickness of 10 mm above or 6 mm underneath the experiment (low energy 
particle shield, LEPS). The top shielding was covered with a blanket of Beta cloth0 and a 
Nomex/ Aluminum/Mylar composite with a thickness of approximately 1.7 mm. By these 
means AMS-01 was protected against mechanical damage and excessive warming from 
exposure to direct sunlight. 

Figure [5TT] (left) shows the definition of the AMS-01 coordinate system. Its origin lies 
in the center of the magnet bore. The z-axis is perpendicular to the tracker and TOF 
planes and points upward, with the Cerenkov counter at negative values of z. The x-axis 
points into the direction of the magnetic field's dominant componentll. 

5.1.1 The Magnet 

A permanent magnet with a nominal bending power of BL^ = 0.14 Tm^ was chosen for 
AMS-01 because of its compact and rigid structure and the low amount of maintenance 
and operating expenses. The magnet had a cylindrical shape with an inner radius of 
558 mm, 95 mm wall thickness and a length of 800 mm along its cylinder axis, resulting 
in a free area of 0.98 m^ and a geometrical acceptance of 0.82 m^sr. It consisted of 6400 
blocks of Nd-Fe-B alloy, sized 2x2x1 inches each, arranged in 64 sectors. The blocks were 
magnetized in such a way that the resulting dipole field was to first order perpendicular 
to the longitudinal axis (z-axis) and parallel to the positive x-axis as shown in Figure 15. 2[ 




Fig. 5.2: x-y cross section through the AMS-01 permanent magnet: arrangement of block 
sectors and AMS coordinate system definition (left); magnetization of blocks and resulting field 
lines (right). 

The magnetic field was not homogenous. As seen in Figure 1531 (left), its dominant 
component reached a fiat maximum of 0.14 T at the center of the magnet and dropped 

^Bcta cloth is a fabric woven from fine quartz threads and impregnated with teflon or sihcone based 
materials |127] . 

^Details about the connection of the AMS-01 coordinate system to the Shuttle and the Earth frames of 
reference are given in [128] . 
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quickly with increasing values of \z\. At the upper and lower edges of the magnet its 
strength decreased to about 64 % of its maximum, while it was down to 14.9 % at z = 
83 cm, the upper end of the detector. Furthermore, was not constant within a cross 
section in the x-y-plane (Figure [531 ^^sht), but grew with increasing values of \x\ up to 
approximately 250 mT. The magnet had very small flux leakage outside its volume and 
a vanishing dipole moment, since the latter, in combination with the geomagnetic field, 
would have implied a torque on the spacecraft. The strength of the magnetic field as a 
function of the space coordinates was measured with high accuracy and is available for 
the offline analysis of AMS-01 data. 




80 60 40 20 -20-40-60-80 [cm] gg eo 40 20 -20-40-60-80 [cm] 



Fig. 5.3: The strength of the x-component Bx of the magnetic field: vertical cross section at 
X = (left) and horizontal cross section at z = (right). White areas denote negative values. 

5.1.2 The Time of Flight Hodoscopes 

The time of flight scintillation counter hodoscopes (TOF) of AMS-01 basically had three 
different tasks: 

• delivering a fast trigger signal to the data acquisition system, 

• measuring the value and direction of an incident particle's velocity, 

• determining the absolute charge of a particle from its energy loss in the scintillators. 

The TOF system consisted of four planes of Bicron BC-408 scintillator paddles, each 
plane covering a total area of roughly 1.6 m^. The planes were grouped in two stations, 
i.e. double layers of scintillator paddles oriented perpendicularly to each other, allowing 
two space points to be measured on the particle trajectory, one above and one below the 
magnet. The distance of two paired planes was 12 cm, the distance from the upper to 
the lower station approximately 130 cm. Thus a particle with (3^1 traverses the time of 
flight system in about 4.3 ns. 

One TOF plane was composed of 14 rectangular scintillator panels, 1 cm thick and 
11 cm wide with a length between 72 and 136 cm to account for the circular shape of the 
plane. In order to avoid dead space, adjacent paddles were mounted with 0.5 cm overlap. 
Each paddle was wrapped in aluminized mylar foil and encased in a 0.6 mm thick carbon 
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Fig. 5.4 '- Schematic view of the upper TOF planes and the LEPS. 

fiber housing. Each scintillator double layer was supported by a 10 cm thick aluminum 
honeycomb panel glued between two aluminum skins of 5 mm thickness |129j . 

A total of six photomultipliers (Hamamatsu R5900), grouped in two triplets, was used 
to read out a single paddle. Both sides of the paddles were terminated with one triplet of 
photomultipliers via 5 cm long trapezoidal light guides. The anode signals from a triplet, 
as well as those from the next to last dynodes, were summed to form two output signals 
for each side of a paddle. 

Prior to the launch of the experiment the time resolution of a single scintillator plane 
was determined from cosmic ray tests |129] to be 126 ps, resulting in an error of 178 ps on 
the total flight time. The position of passage of a particle through the scintillators can 
be determined with an accuracy of 1.8 cm along a panel. 

The charge measurement is performed using the time information, based on a time- 
over-threshold method. Though not optimized for this task, the TOF system was thus 
capable of distinguishing particles with charge \Z\ = e, \Z\ = 2e and \Z\ > 2e with a 
purity better than 1 % [130] . 

5.1.3 The Silicon Strip Tracker 

For the measurement of position and energy loss of particles AMS-01 featured a high 
quality tracking device based on silicon microstrip sensor technology |131j . It was designed 
to provide a position resolution of 10 /im in the bending plane and 30 fim in the non- 
bending plane of the magnetic field. 

Tracker Layout and Electronics 

The layout of the tracker is displayed in Figure 1531 (left). The silicon was arranged in 
6 planes roughly 20 cm apart (see Table [5?T]) . Four of the planes were located inside the 
magnet bore, enclosed in a carbon fiber cylindrical shell. Each plane was supported by a 
12 mm thick aluminum honeycomb disk dressed with 220 fim carbon fiber skins. Above 
and below the magnet, two additional layers were mounted on stiffer disks (4 cm thick 
with 0.7 mm cover) for mechanical stability. 



Layer no. 


1 


2 


3 


4 5 


6 


z position [cm] 


50.98 


29.19 


7.79 


-7.82 -29.18 


-50.98 



Tab. 5. 1 : The z position of the tracker layers in the AMS coordinate system. 
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Fig. 5.5: Schematic view of the siUcon strip tracker and its support structure (left); picture 
of the fifth tracker layer during detector disassembly (right) . 

The basic tracker elements were the 300 fim thick silicon microstrip sensors of size 
40.14 X 72.04 mm^. The sensors were double sided with parallel strips implanted in the 
n+ bulk silicon on both sides. Strips on opposite sensor sides were perpendicular to each 
other [132j . On the p-(junction-)side, which measures the y-coordinate, the implants had 
a pitch of 27.5 /xm with every second implant covered with an aluminum readout strip on 
the surface. On the ohmic n-side the implantation pitch was 26 fim with every implant 
aluminized. Every fourth implant was read out on both sides, resulting in a readout pitch 
of 110 /im for the p-side and 208 fim for the n-side. 

Between 7 and 15 sensors were arranged in linear structures, the so called ladders, of 
different length to account for the overall circular shape of the layers. Within a ladder, the 
p-side readout strips of all sensors were daisy chained with wirebonds and connected to 
the p-side frontend hybrid. The n-side strips were connected to the n-side frontend hybrid 
through a metalized kapton foil, directly glued to the sensors. The kapton foil incorpo- 
rated 2 groups of 192 lines each, with adjacent sensors connected to a different group. As 
a consequence of this bonding scheme |133] (see Figure 15. 6p and depending on the number 
of sensors in a ladder, between 4 and 8 strips on different sensors were connected to the 
same readout channel, thus leading to an ambiguity in the assignment of tracker hits 
to the n-side strips. The two hybrids were attached back-to-back at each ladder's very 
end, close to the magnet wall, and were perpendicular to the sensor planes. They were 
equipped with 64-channel VAhdr readout chips |134] . which performed charge-sensitive 
signal preamplification and shaping and also provided a sample-and-hold stage, buffer 
and multiplexer for sequential analog data output. In the onboard readout electronics, 
the output was digitized by fast 12-bit ADC after amplification. Further downstream, 
digital signal processors performed pedestal subtraction, noise determination and clus- 
tered the strip signals, so that compressed hit information was transmitted to the data 
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X 



Fig. 5.6: The principle of n-side strip connection to the frontend electronics as a cause for 
ambiguities in the hit reconstruction. Connection via one kapton track group only is displayed. 



acquisition of the experiment (zero suppression). For clustering, strips with signal am- 
plitudes more than three times above their noise level cxped were used as seed strips and 
grouped with neighboring strips having amplitudes above 1 aped to form clusters of at 
most 5 strips |131j . The cluster position was determined from the center of gravity of the 
strip signal amplitudes. 

As evident from Figure 15.51 (right), the tracker planes were not fully instrumented. 
From the full number of 1912 sensors, 778 were finally assembled |135j . covering a total 
area of more than 2m^. In this configuration, the tracker had an acceptance of 0.31 m^sr 
for particles traversing four or more layers of instrumentation [132] . 

Tracker Characteristics and Performance 

Using the framework of the offline software, the characteristics and the performance of the 
tracker were determined with single-track events from flight and from test beam data. Two 
different methods were used for track fltting, based on an iterative Kalman fllter approach 
or on the numerical integration of the equation of motion of particles in a magnetic fleld. 
For singly charged particles. Figure 15.71 [left) shows the rigidity resolution AR/R as a 
function of rigidity, for tracks with 4, 5 or 6 measured points from flight and CERN-PS 
test beam data, in comparison with Monte Carlo results [136] . Multiple scattering limits 
the rigidity resolution to a level of 9 % in the low rigidity range up to approximately 
10 GV. For higher energies, the resolution deteriorates rapidly as a consequence of the 
finite magnetic field and the sensor spatial resolution. 

On the right side of Figure 15.71 the average energy loss per tracker plane is displayed 
for protons and helium nuclei from fiight data |131] . The energy loss is derived solely 
from p-side strip amplitudes and corrected for path length in the silicon. Fits of a sum of 
Landau and Gauss functions to the distributions give a most probable energy deposition 
of 96keV for protons and 397 keV for helium nuclei. 

The detection efficiency for a single plane is determined by the number of bad channels 
in the silicon, as well as by the fact that adjacent ladders were mounted on the support 
structures without overlap. The average total efficiency of the individual tracker planes 
varied between 80 % and 85 %, depending on the signal to noise ratio, which is a function 
of tracker temperature and thus of time. The contribution from bad channels - which 
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Fig. 5.7: The tracker rigidity resolution as a function of rigidity |136j determined for well 
defined tracks from flight and test beam data and compared with Monte Carlo simulations 
(left); truncated mean of energy loss in the tracker planes |131j for protons and helium nuclei 
(hatched histogram) from data, based on p-side strip measurements (right). 

may be due to high noise level, defective readout electronics, or non-linear signal response 
- sums up to approximately 12 % of the inefficiency. 

During the Space Shuttle flight there were deviations of the tracker support structure 
position from its nominal value in the order of 30 fim due to temperature variations. These 
deformations were monitored with an onboard infrared laser system |137j . The tracker 
was flnally aligned through post-flight metrology and the analysis of high-rigidity tracks 
from flight data |138j . 

5.1.4 The Aerogel Threshold Cerenkov Detector 

Mounted underneath the lower time of flight scintillator layers, the Aerogel Threshold 
Cerenkov Detector (ATC) provided separation of antiprotons from electrons or protons 
from positrons for momenta below 3.5 GeV/c. 




Fig. 5.8 : Cross section of the Aerogel Threshold Cerenkov Detector and detail of the structure 
of an aerogel cell |139j . 
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The ATC was composed of aerogel cells with a volume of llx 11 x8.8 cm^, as shown in 
Figure EiHl Each cell was filled with 8 slabs of 1.1 cm thick aerogel blocks and surrounded 
by 3 layers of teflon, each 250 /im thick, to reflect Cerenkov photons back into the cell 
volume. At the bottom end, a photomultiplier was attached to the cell via a light guide. 
The efficiency of the photomultiplier reached a maximum at wavelengths of approximately 
420 nm. In order to reduce photon losses due to scattering and absorption, a 25 fim 
wavelength shifting layer was inserted between the aerogel slabs in the middle of each 
cell. 

The 168 aerogel cells were grouped into two layers of 8x10 and 8x11 cells respectively, 
and were enclosed in a thin carbon fibre structure. Bolted to a 5 cm thick honeycomb 
plate which was glued into an aluminum frame, the layers were attached to the support 
structure on which the AMS-01 experiment was mounted in the Space Shuttle payload 
bay. 

The data from the ATC are not used in this analysis, since they provide no background 
suppression in the high momentum range above 3.5 GeV/c. A detailed description of the 
device is given in [139] . 

5.1.5 The Veto Counters 

The 16 veto or anti-coincidence counters (ACC) were mounted on the inner magnet sur- 
face, thus surrounding the tracker. Their purpose was to provide detection of particles 
which entered or exited the detector volume through the magnet wall, thus typically indi- 
cating multi-particle events involving 5-ray creation or the presence of low quality tracks. 
Each counter consisted of a 1 cm thick scintillator paddle 80 cm high and 20 cm wide with 
implanted wavelength shifting fibers. At both ends of the paddles, the fibers were viewed 
by a single photomultiplier. 

5.1.6 Material Budget 

Depending on direction of flight, particles had to traverse a certain amount of material 
before they could reach the tracker. This fact has important implications for this data 
analysis. Additional material induces multiple scattering of particles and hence reduces 
the spatial and angular resolution of reconstruction algorithms. On the other hand, it 
increases the probability of processes such as bremsstrahlung and gamma conversion, 
which are the focus of this work, and therefore improves the reconstruction efficiency. 

Table 15.21 summarizes the dominant contributions to the overall material budget from 
particular components and subdetectors of AMS-01 outside the tracker. While the thermal 
shielding and the LEPS were located on top of the experiment, the aerogel was in the 
lowermost subdetector of AMS-01 and was located below the tracker. As a consequence, 
for particles entering AMS-01 from the top (downward), the total amount of material to be 
traversed before reaching the tracker was equivalent to 18.2 % of a radiation length, while 
for particles entering from the bottom (upward) it amounted to 19.2%. Each tracker 
layer, including the silicon and the support structure, contributed another 0.65 % of a 
radiation length to the material budget |131j . 
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Component 




Ref. 


Thermal shielding 


3.2 


\m\ 


LEPS 


3.9 


|140| 


TOF double layer 


9.7 


m 


TOF double layer support 


1.4 




Aerogel 


8.0 


|141| 



Tab. 5.2 : Contribution to the overall material budget from the components of AMS-01 outside 
the tracker, in percent of a radiation length. 



In addition to the components of the experiment itself, additional objects have to be 
considered for the material budget for upward going particles, namely the Space Shuttle 
and the AMS-01 support structure in it. As shown in Figure the mounting frame of 
AMS-01 (Unique Support Structure, USS) consisted of aluminum box beams and truss 
elements, some of which were within the acceptance of the detector. The 4 mm thick 
outer aluminum wall was covered with low density silica fiber tiles |142j . forming the heat 
shield. Besides, the floor contained trusses, as well as hydraulic and electric circuitry. 
Additionally, the inner surface of the payload bay was lined with a protective blanket 
made of composite fabric material. Precise information about the structure of the Space 
Shuttle's payload bay floor, which would allow an exact determination of the material 
distribution underneath AMS-01, is not available. The amount of this additional material 
is estimated to be 4.5 % of a radiation length from the specifications of dimensions and 
weight of the Space Shuttle fuselage |142] . 

5.2 Trigger and DAQ Livetime 

The trigger chain of AMS-01 was divided into three successive stages: the hardware- 
based Fast Trigger, and the Level-1 and Level-3 stages, which were implemented in the 
data acquisition Each trigger stage definitely rejected an event if its particular 

requirements were not fulfilled |143j . For trigger efficiency studies, a small fraction of 
events - approximately 1 in 1000 - was accepted after a positive decision of the Fast 
Trigger only (prescaled events). 

The Fast Trigger signal was derived solely from the photomultiplier triplets (summed 
anode signals) viewing the TOF scintillators. A Fast Trigger was generated if at least 
one scintillator panel signaled a hit from one or both end triplets in each of the four TOF 
planes. The typical time for the Fast Trigger stage to generate the signal was 100 ns. If 
after another 200 ns none of the following trigger requirements were fulfilled, the event 
was rejected. 

The first trigger condition of the Level-1 stage was called matrix condition and ac- 
counted for the fact that the tracker layers were not fully equipped. By comparison with a 
correlation matrix, those combinations of TOF scintillator panels on layers 1 and 4 which 

•^A Level-2 stage was not implemented. 
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were incompatible with a particle traversing the tracker instrumentation were identified 
and the corresponding events were rejected. The second trigger condition of the first stage 
demanded the absence of signals in any of the veto counters. 

The Level-3 stage assured the coincidence of the two photomultiplier triplets on either 
side of the TOF panels used for the fast trigger decision. Furthermore, it performed a 
clustering of the TOF hits, forming clusters either from one hit in a single panel or from 
two hits in two adjacent panels. A single cluster was then required in each of the first 
and the last TOF planes in coincidence with at least one additional cluster in the second 
and/or the third plane. These TOF clusters defined a 6 cm wide corridor in the bending 
plane. The Level-3 stage retained the event, if at least three tracker hits were found within 
this corridor with a signal to noise ratio above 4 [136] . Events with multiple clusters on a 
single TOF plane were rejected |144j . The latter trigger criterion has a negative impact 
on this data analysis, since it deteriorates the acceptance for the signal process towards 
lower particle momentum (see § I7.2.ip . 

An additional feature of the Level-3 stage was active during the first seven hours of 
data taking |144j . Its main purpose was the rejection of low energy particles with \Z\ < 3. 
The algorithm relied on the measurement of the average residual of tracker clusters with 
respect to a straight line joining the uppermost and the lowermost cluster. Events were 
accepted if the average residual was below 300 /xm. 

The total dead time per trigger operation amounted to 85 /xs [IB]. In case an event was 
rejected by the Level-1 stage, it took 7/is until the data acquisition system (DAQ) was 
again operational. Consequently, the fraction of time during which the DAQ was active 
[DAQ livetime) significantly depended on the particle rate, which itself varied with the 
position of the detector in the geomagnetic coordinate system. At low magnetic latitude, 
particle rates were small and the readout rate reached up to L6kHz. Except over the 
South Atlantic, where the geomagnetic field is exceptionally weal0 resulting in very low 
readout rates due to a large amount of trigger dead time, the minimum readout rate was 
100 Hz. 

5.3 AMS-01 on the Space Shuttle 

The AMS-01 experiment was fiown on the Space Shuttle Discovery in the course of the 
STS-91 missio ifl from June 2nd to 12th, 1998. Fi gure 15.91 shows the mounting of AMS-01 
in the spacecraft's cargo bay. During its fiight the Discovery orbited the Earth at altitudes 
between 320 and 390 km with an orbital period of 93 min and an inclination of ± 51.7 de- 
grees. From June 4th to 8th, for a total of 95 hours, the Space Shuttle was docked to 
the Mir space station. During docking, parts of the Mir, such as the Soyuz-TM and the 
Priroda modules, were within the acceptance solid angle of the AMS-01 experiment |145j 
(see Figure Em]) • 

The 184 hours of data taking were divided into several periods. As seen in Figure [3.101 
the pointing angle of AMS-01 - the angle of the positive z-axis with respect to the zenith 
direction - varied with time. Before docking to the Mir, the pointing angle was kept 

phenomenon known as the South Atlantic Anomaly (SAA) 
^STS: Space Transportation System 
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constant at 45 degrees. After decoupling from the space station, the experiment was 
operated at consecutive angles of 0, 20, 30, 45 and 180 degrees. During the latter phase, 
the z-axis of AMS-01 pointed towards the Earth for the measurement of outbound particles 
emitted from the atmosphere (albedo phase). While the Discovery was docked to the Mir, 
the pointing angle could not be kept constant; it oscillated between roughly 40 and 140 
degrees with the Mir orbital frequency of about 0.65 h~^. Several short loss-of-signal 
periods occured during the flight, partly because the electronics had to be shut down due 
to overheating when the shuttle was pointing towards the zenith |146j . 



The AWS Experimerrt 
Coortiinote Sy-stcm: 
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Fig. 5.9: The mounting of AMS-01 in the Space Shuttle's cargo bay (left); photograph of the 
Space Shuttle taken from the Mir space station (right). 
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Fig. 5.10: The pointing angle of AMS-01 as a function of time. 
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Fig. 5.11: Schematic view of the Space Shuttle docked to the Mir space station |147| . The 
Soyuz-TM and the Priroda modules were within the acceptance of the AMS-01 experiment. 
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6 Positron Identification through 
Bremsstrahlung Conversion 



6.1 The Challenge of Positron Measurements 

The main challenge of high energy cosmic ray positron measurements is the suppression 
of the vast proton background. As already stated in chapter HI the flux of cosmic ray pro- 
tons exceeds that of positrons by a factor of 10^ in the momentum range of 1-50 GeV/c. 
Besides the fact that both particle species have the same charge, their energy deposition 
in traversed material becomes comparable at energies beyond a few GeV. Thus, the iden- 
tification of positrons must rely either on techniques based on measurement of the Lorentz 
factor, as realized for example in Cerenkov counters or transition radiation detectors, or 
on high-resolution calorimetry. 

However, as described in the last chapter, AMS-01 was not equipped with a calorimeter 
and its ATC subdetector provided a sufficient proton rejection only for momenta below 
3.5 GeV/c. In order to extend the energy range accessible to the experiment, a completely 
different approach has been chosen. It relies on the identification of bremsstrahlung 
emission through photon conversion, the latter also referred to as pair production. 



6.2 Conversion of Bremsstrahlung Photons 



Both bremsstrahlung and photon conversion cannot occur in vacuum, but require the 
presence of a Coulomb field, predominantly that of a heavy object like a nucleus, so that 
momentum and energy conservation are satisfied. They are closely related electromagnetic 
processes whose energy-angle distributions can be calculated with the Bethe-Heitler for- 
malism. For bremsstrahlung from electrons or positrons and in the limit of high energies 
this distribution is given by the cross section |148] 
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where Eq and E denote the energy of the lepton before and after the radiation process, 
k = Eq — E the energy of the photon, Z the charge number of the scattering nucleus, 
m the electron mass and c the speed of light. C is a constant of the order of 10^, which 
expresses the screening effect of the atomic electrons bound to the nucleus. The variable 
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X is equal to the reduced angle 76*, 7 being the lepton's Lorentz factor and 6 the angle of 
the radiated photon with respect to the incident lepton. 

The bremsstrahlung matrix element is related to those of pair production by the sub- 
stitutions k ^ —k and p ^ —p, where p is the four-momentum of the incident particle 
in bremsstrahlung emission or the four-momentum of one of the pair particles in pair 
production. Thus, the energy-angle distribution for pair production has a mathematical 
form similar to that for bremsstrahlung, with 6 being the angle of one of the leptons with 
respect to the incident photon |149j . 

The above cross section (16.11) depends inversely on the square of the incident parti- 
cle mass. Hence it follows that bremsstrahlung emission is suppressed by a factor of 
more than 3 ■ 10^ for protons with respect to positrons. This important property of the 
bremsstrahlung process has a very beneficial implication: identification of the radiated 
photon, in conjunction with a measurement of the charge sign of the incident radiating 
particle, determines the latter almost certainly as a positron (or electron). This principle 
is the basis for the proton background suppression as performed in this analysis. 

6.2.1 Event signature 
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Fig. 6. 1 : Schematic view of a converted bremsstrahlung event caused by a downward going 
positron. 

Figure 16.11 shows the topology of an event with a converted bremsstrahlung pho- 
ton. Here, a primary positron enters the detector volume from above and emits a 
bremsstrahlung photon in the first TOF scintillator layer. The photon then converts 
into a secondary electron-positron pair in the second TOF layer in this example. If the 
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conversion process takes place above or in the upper part of the tracker, the tracks of 
three particles can in principle be reconstructed. 

Eq. (16.1 p can be integrated over x to obtain the root-mean-square angle 6rms of the pho- 
ton with respect to the incident lepton, yielding 9rms = (i{Eq, k, Z)/7-ln7 for high energies 
[150] . Since q{EQ,k,Z) depends primarily on k/Eo and is always of the order of unity, 
it follows that Orms ~ 1/7. As a consequence, in the GeV energy range, bremsstrahlung 
photons are largely emitted under angles of 10""^ rad or smaller. The same applies to the 
opening angles of the electron positron pairs from pair production. In a tracking detector, 
these values fall below the typical limit of track reconstruction accuracy given by multiple 
scattering and thus are practically equal to zero. 

Furthermore, it can be seen from eq. fl6.1l) that the spectrum of the emitted photons is 
soft, since the cross section decreases with 1/k. Because of the low fraction of momentum 
which is typically carried away by the photon, the secondary particles have lower momenta 
than the primary. Therefore, in the bending plane projection of the spectrometer, the 
secondaries tend to form the left and right tracks, while the primary remains in the middle 
(see Figure 16. ip . It must be pointed out that in the GeV energy range the fraction of 
momentum given to each of the secondary particles is approximately uniformly distributed 
between and E^. In the limit of high energies, the distribution of the ratio of the two 
secondary momenta becomes slightly asymmetric, with electron-positron pairs having 
considerably different momenta being favored. 

6.2.2 Background 

The background to positron identification through bremsstrahlung conversion is caused 
by protons undergoing hadronic reactions in the material of the experiment, as well as by 
electrons with wrongly reconstructed momentum sign. 

Mesons produced in hadronic reactions involving protons can mimic the three-track 
signature of converted bremsstrahlung events. These mesons are almost exclusively pions; 
the cross section for the production of kaons or other mesons is in comparison more than 
one order of magnitude smaller [151j . Several processes contribute to this background 
and can be grouped in two categories, where p, n denote an incident proton, or a proton 
or neutron in the scattering nucleus (cross sections are taken from |152] ): 

• processes involving charged pion production, dominantly pp — ppir'^ir' and pn 
pnvr+Tr", with cross sections of approximately 2mb above a threshold of 3 GeV 
incident proton energy, and 

• processes involving neutral pion production, such as pp — > ppn'^ or pn — > pmr^, 
with cross sections of approximately 2mb or 4mb, respectively, above a threshold 
of 1 GeV incident proton energy. 

In case of neutral pion production, the vr'^'s decay almost exclusively into photon pairs 
after a mean flight path of cr = 25 nm. If one of the photon remains undetected, while 
the other one converts into an electron positron pair, the resulting three-track event 
signature is that of a converted bremsstrahlung event. Due to the procedure for vertex 
reconstruction, discussed in § 16.3.41 the fact that - in contrast to the signal process - 
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only one vertex is present in the above hadronic interaction cannot be used to effectively 
discriminate against the hadronic background. 

The material of the experiment, in which the relevant reactions take place, mainly 
consists of carbon and hydrogen atoms. According to the optical model |153] , the proton- 
proton cross sections cTpp can be applied to proton-nucleus scattering when appropriately 
scaled as (jpN = Upp^", with A being the atomic weight of the target nucleus, and a = 
0.77 for carbon. Hence, for a scintillator of the TOF system with a thickness of 1 cm, 
the probability for a proton to undergo one of the above reactions is of the order of 
1 %o for both charged and neutral pion production. However, the probability for such 
a reaction to give rise to a background event is much smaller due to several reasons. 
The spectral distributions of the secondary pions show distinct peaks at 100 MeV and 
decrease steeply towards higher energies [154] . The resulting momenta are typically very 
low and tend to prevent a hadronic event from exhibiting a three-track signature due to 
the strong deflection in the magnetic field. On the other hand, in contrast to converted 
bremsstrahlung events where the secondary tracks emerge from the vertex at small angles 
and are thus ordered according to their charge sign by the spectrometer magnetic field, 
secondary particles from hadronic interactions are emitted under large angles with respect 
to each other and may be any one of the three tracks with almost equal probability. 
Background from processes involving neutral pion production is further suppressed by the 
requirement that exactly one of the two photons from pion decay must convert. Let p ~ 
7/9Xo denote the pair production probability in a material with a thickness corresponding 
to a fraction Xq of a radiation length |155j . Then the suppression factor is {p{l —p)}~^ = 
8.33 for Xq = 0.18, in the case of downward particles. 

The average transverse momentum (pi) of mesons stemming from proton-proton in- 
teractions is (pt) ~ 0.35GeV/c, rather independent of the collision energy |156] . As 
a consequence, the angles 6 between particles emerging from interactions in the GeV 
energy range are large compared to those in bremsstrahlung conversion. For example, 
(6) ^ 33 • 10~^rad at 10 GeV. This important characteristic of hadronic interactions al- 
lows to effectively discriminate against the proton background. 

Determining from first principles the background from electrons with wrongly recon- 
structed charge sign is naturally harder due to the rather complex behavior of track 
reconstruction algorithms. As Figure 15.71 on page [59] shows, the rigidity resolution of the 
tracker has a constant value of about 9 % in the range 1-9 GV and then rises linearly 
with increasing rigidity. For high energy particles, the resolution rapidly worsens so that 
the probability of bad reconstruction increases. If the resolution is approximated by the 
width 0" of a normalized Gaussian distribution G{R, R, a) of the rigidity R around a true 
mean value R (which is negative for electrons), then the probability Pw of wrong charge 
sign reconstruction is given by the integral 



Due to the quickly falling edges of the Gaussian distribution, the value of is practically 
zero for low rigidities. At 50 GV it amounts to approximately 1 %o and hence the amount 
of misidentified electrons remains small compared to the positron to electron ratio. 
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Apart from these rather coarse estimations, the properties of the background sources, 
as well as their individual contributions to the signal event sample, must be precisely 
determined by Monte Carlo simulation. In § 16.4.21 it will be shown that in the course of 
the analysis the probability for a proton to create a background event can be reduced to 
less than 10~^, and the contamination of the positron sample from badly reconstructed 
electrons remains below the percent level. 

6.3 Event Reconstruction 

The task of the event reconstruction process is to identify candidates for converted 
bremsstrahlung events and reconstruct the event geometry and particle momenta. In 
order to obtain the highest possible selection efficiency, it is mandatory to apply sophis- 
ticated track and vertex finding algorithms which are particularly customized for the 
signature of converted bremsstrahlung. In a first step, a fairly permissive preselection 
algorithm extracts signal candidates from the AMS-01 data set. Two independent track 
finding algorithms, each of them optimized for a particular event topology, are separately 
applied to the preselected event sample. In case that three particle tracks are found, their 
corresponding momenta are determined, and the bremsstrahlung and conversion vertices 
are spatially reconstructed. As argued in § 16.2.11 the middle track is considered to stem 
from the primary particle, while the left and right tracks represent the secondary electron 
positron pair. To account for the asymmetric geometry of the detector along its z-axis, 
event reconstruction and analysis are performed separately for particles traversing the 
detector from top to bottom (downward) and from bottom to top (upward). 

Before the analysis, the raw data recorded by AMS-01 were preprocessed |157j . This 
procedure included, amongst others: 

• reconstruction of the spatial coordinates of the tracker hits and their errors, 

• reconstruction of the energy deposits and relative transit times of particles in the 
TOF scintillator layers, and 

• calculation of the trigger livetime, i.e. the fraction of time over which the trigger 
electronics of the experiment were accepting signals (see § I7.2.2p . 

6.3.1 Preselection 

The event reconstruction algorithms developed for this analysis are rather complex and 
require substantial amounts of computing time because of their combinatorial character. 
The task of the preselection is to considerably reduce the size of the AMS-01 data, so 
that the algorithms can operate on a reasonably small sample of event candidates. 

In a first step, all hits in the silicon strip detectors of the tracker are projected into 
the bending plane for one dimensional clustering. On a particular tracker layer, a cluster 
seed is defined by a single arbitrary tracker hit. Subsequently, further hits are added to 
an existing cluster if their distance from the barycenteio of this cluster along the y-axis is 

^The barycenter is the weighted mean of the hit coordinates, with the hit amphtudes as weights. 
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smaller than 100 fim. If a hit cannot be added to any existing cluster, it is used to start 
a new cluster. In this procedure, the x-coordinates of the tracker hits are not regarded, 
and the clusters are defined in the bending plane (y-z plane) only (see Figure [U75|l . Thus 
the hits in a given cluster may spread over the whole length of the particular tracker layer 
along the x-axis. Most of them are virtual hits caused by the tracker bonding scheme (see 
§ I5.1.3P . This ambiguity is resolved by the track finding algorithms in a subsequent step 
(see § 16.3.21) . Moreover, there may be clusters comprised of a single hit only. 

For further analysis, a minimum of 8 tracker clusters are required. Events are selected 
in which at least two of the six layers of silicon detectors contain exactly three clusters 
each (triplets). By these means, a reduction of the data size by a factor of 25 is achieved. 



6.3.2 Track Finding 

AMS-Ol was designed for single particle detection, thus the software package for data 
analysis provided by the AMS-Ol collaboration is optimized for the reconstruction of 
single-track events. The track finding algorithms implemented in this package have a 
small reconstruction efficiency for multi-track events, resulting in low event statistics when 
applied to a rare process such as bremsstrahlung conversion. As a consequence, alternative 
reconstruction methods had to be developed for this analysis. Two different track finding 
algorithms - in the following referred to as the successive and the combinatorial algorithm 
- are independently applied each to a subsample of the preselected events. Their purpose 
is to establish if an event exhibits a three-track signature and, if so, find the correct 
assignment of the tracker clusters to the three tracks. Each algorithm is optimized for a 
particular cluster topology and returns a set of parameterized track information. 

With t denoting time, and p and v being the momentum and velocity of a particle of 
charge e, its equation of motion in a magnetic field B is 

^ = e(vxB). (6.3) 

In a homogeneous field, the solution of this equation is a helix trajectory. Since the 
magnetic field within AMS-Ol is inhomogeneous, eq. (16. 3p is preferably solved numerically. 
Nevertheless, the resulting trajectories can still be well approximated by a helix. The 
latter is fully parameterized with a set h of five numbers, defined on any given plane 

z = zq as 

dx dy 1 



dz dz p 



(6.4) 



An iterative fitting algorithm [158j . particularly adapted for the AMS-Ol experiment [159] . 
is used to determine h from a set of point measurements i of a particle's trajectory. The 
algorithm is based on minimization of the quantity 
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with Xi,yi being helix parameters on the given z plane and Xi,yi the measured hit coor- 
dinates on that plane with the corresponding measurement errors (Tx,ay. The algorithm 
returns h at the z-position of the first measurement in the direction of flight. The mini- 
mum rigidity that can be reconstructed is 100 MV. 

The successive track finding algorithm 

The successive track finding algorithm is customized for events in which at least three 
tracker layers are found with exactly three clusters on each of them. Its working principle 
is illustrated in Figure 16. 2[ Since particle tracks diverge in the magnetic field, the triplets 
- clusters on layers with exactly three clusters - are required to have increasing cluster 
distances along the z-axis in the direction of flight. With three particles having traversed 
the tracker, the clusters in the triplets can directly be assigned to a left, a middle and a 
right track segment (initial track seeds). 

Subsequently, starting from the uppermost layer, further clusters on the other layers 
are tentatively added to each of the tracks. A cluster is finally appended to the track 
which yields the lowest according to eq. fl6.5l) . if the value of is below a threshold 
of 30 (the fitting procedure is described in the next paragraph). In this case, the helix 
parameters of this track, including the new cluster, are updated. As the tracks grow, as 
many clusters as possible are attached to them. Clusters caused by noise or the occurrence 
of 5-ray electrons are mostly incompatible with any of the tracks and remain isolated. 

For the fitting procedure mentioned above, it is necessary to resolve the ambiguities 
in the x-coordinates of the clusters, caused by the clustering only in the bending plane 
projection. For this, a linear regression procedure |160j is carried out to obtain a straight 
line through the hits in the TOF system (see Figure [HIS]) • In the non-bending plane (x-z- 
plane), a corridor of 10 cm width is defined around the straight line and only hits within 




Fig. 6.2: Principle of the successive track finding algorithm in the bending plane projection. 
Initial track seeds are created from seed clusters on layers with exactly three clusters. Further 
clusters are successively added to the tracks, while noise clusters, which do not fit to any track, 
remain isolated. Missing clusters occur in parts of the tracker which are inefficient or not 
instrumented. 
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this corridor are retained. All possible combinations of remaining hits in the track - with 
at most one hit from each cluster - are then tentatively fitted, and the helix parameters are 
updated from the combination which yields the lowest value of according to eq. fl6.5l) . 
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Fig. 6.3: Display of a three-track event, projected in the bending plane (left) and in the 
non-bending plane (right). Only the TOF system and the tracker are shown. The height of 
the hit boxes is proportional to the hit amplitudes, while their width represents the position 
measurement errors. On the right side, the hit selection corridor is plotted. 



The combinatorial track finding algorithm 

In order to increase the efficiency for converted bremsstrahlung events, a generalized 
algorithm has been developed for the treatment of events that feature only two tracker 
layers with exactly three clusters each |161j . It is based on a combinatorial approach to the 
track finding problem. The basic idea is to examine every possibility (track hypothesis) of 
arranging the tracker clusters into three tracks and determine those combinations yielding 
the lowest value of according to eq. fl6.5p . These track hypotheses have to pass several 
tests in order to reject wrongly assigned tracks. Additionally, two tracks may share one 
single cluster. The algorithm leads to an improvement of the efficiency by 40 %. A detailed 
description of this method is given in [161] . 



6.3.3 Vertex Reconstruction 

Vertex reconstruction is based on parallel backtracing of tracks through the magnetic field 
of the detector using the functionality of the GEANT3 package (routine GRKUTA) [162] . 
For a given track of a particle i with charge q and momentum p, a virtual particle i with 
charge —q and momentum — p is started at the first tracker cluster found in the fiight 
direction of particle i. The propagation of i is then performed by numerically solving the 
equation of motion in steps of 100 /im, so that i will follow the trajectory of the original 
particle i backward against its original flight direction. 
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The vertex finding for tlie case of downward going particles is performed as follows. 
The left and right tracks may have their starting cluster on different tracker layers. As a 
first step, the track starting at the lower layer is backtraced upward until it reaches the 
starting layer of the other track. Both tracks are then backtraced in parallel until their 
projected distance along the y-axis reaches a minimum. The vertex is then defined as the 
barycenter of the track points at the z-coordinate of closest approach of the tracks. In 
case the tracks intersect in the bending plane projection, the intersection point is taken 
as the vertex. The vertex finding for upward going particles is exactly similar, but with 
all directions reversed in z. 

Finally, the four-vector of the photon is reconstructed from the four-vectors of the left 
and right tracks at the vertex position. Using the algorithm described above, the vertex 
of the photon and the middle track is then computed. In this case, the photon is simply 
propagated along a straight line. 

6.3.4 Reconstruction Quality and Monte Carlo 

The quality of the reconstruction algorithms is verified with 16.8-10^ electron and positron 
events from a complete Monte Carlo simulation of the experiment using GEANT3. The 
incident momentum spectrum of these Monte Carlo particles is proportional to p^^ and 
ranges from 1 to 60 GeV/c. For background studies an additional 170 ■ 10^ proton Monte 
Carlo events have been generated in the momentum range 6-150 GeV/c. All Monte Carlo 
particles are injected both from above and below into the acceptance solid angle of the 
experiment with a uniform spatial and angular distribution. 




Fig. 6.4 '- Comparison of generated and reconstructed momentum (sum of three track mo- 
menta) for downward (left) and upward (right) going positrons. 
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Fig. 6.5: Comparison of generated and reconstructed bremsstrahlung photon momentum 
(left) and the absolute error in its direction (right) for downward and upward going particles 
combined. 

Figure 16.41 shows the ratio of generated and reconstructed momentum of converted 
bremsstrahlung events as the sum of the three track momenta for downward (leit) and 
upward (right) going positrons. A Gaussian fit results in a momentum resolution of 
approximately 14 % and 15 % for the downward and upward case, respectively. This 
resolution is quite comparable to that for single-track events in the energy range of 10 GeV 
and above |131l 1163] , where the reconstruction algorithms have their peak sensitivity (see 

§[72]). 

The tails towards higher relative momenta are caused by additional bremsstrahlung 
photons which are emitted from the primary before it enters the tracking volume. If such 
a photon remains undetected, the measured momentum of the primary will be lower than 
the true initial momentum. This effect is especially apparent in the case of upward going 
particles due to the high amount of material they have to traverse before reaching the 
tracker. 

As displayed in Figure 16. 5^ the properties of the bremsstrahlung photon can be parti- 
cularly well reconstructed. A Gaussian fit to the photon's relative momentum results in 
a momentum resolution better than 9 %, while the absolute error in its direction has a 
standard deviation below 9 mrad. The secondary tracks from which the photon is recon- 
structed typically have lower momenta than the primary particle, resulting in a better 
momentum resolution. Additionally, in the case of the photon, the direction informa- 
tion of both secondary tracks enters into the reconstruction, thus increasing the resulting 
precision. 

The spatial reconstruction of the vertices is rather poor due to the low angles under 
which the particles emerge from them. Figure [6^6] (leR, middle) shows the absolute error of 
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Figf. 6.6: Absolute vertex position reconstruction error zmc — -Zrcc along the z-coordinate (left) 
and in the x-y-plane (middle), schematic view of the vertex finding in case of intersecting tracks 
(right). 



the reconstructed vertex positions, calculated from downward going Monte Carlo electrons 
and positrons. For the conversion vertex, the resolution of the z-position is not better than 
10 cm. Furthermore, a mean displacement of 8 cm is observed, i.e. the vertex is typically 
reconstructed too close to the center of the detector. The reason for this displacement 
lies in the geometrical characteristics of the vertex finding algorithm. As displayed in 
Figure [nS] (I'ight), small errors in the track directions may lead to intersecting tracks in 
the bending plane projection, always resulting in a vertex position which is shifted in 
the flight direction of the particle. For the same reasons, the bremsstrahlung vertex is 
reproduced with an accuracy of 16 cm and a displacement of about 13 cm in the flight 
direction. For the conversion vertex, the width of the error distribution in the x-y-plane 
(Figure 16.61 middle) is mostly a consequence of the z-errors when tracks are tilted with 
respect to the AMS-01 z-axis. Since the opening angles of the tracks vary only slightly 
with vertex position, the low quality of the vertex position reconstruction has almost no 
implications for the analysis. 



6.4 Analysis 

Analysis and suppression of background mainly rely on the evaluation of the topology and 
of the geometrical properties of the reconstructed events and are therefore based on data 
from the tracker. Additionally, cuts are applied on data from the TOF system. However, 
substantial parts of the analysis deal with measures to account for the environmental 
circumstances under which the AMS-01 experiment was operated, especially the effect of 
the geomagnetic field. 
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-Fig. 6.7: Flight time along the z-axis measured with the TOF system for downward (top) 
and upward (bottom) going electron and positron candidates from data and Monte Carlo. The 
Monte Carlo distributions are scaled to the data by integral. 



6.4.1 Basic Cuts 

Several simple cuts are applied to the data in order to eliminate wrongly reconstructed 
events and a large portion of the background: 

• As mentioned in § 16.3.21 the reconstruction of tracks with momenta lower than 
100 MeV is not reliable. Tracks reconstructed with momenta below this value are 
tagged by the track fit algorithm and events containing them are rejected. 

• Due to the deflection in the magnetic field and in conjunction with the small emission 
angles, the charge signs of the secondaries are exactly constrained depending on the 
direction of incidence. For downward going particles the left track is required to 
have positive charge sign, while the right track must be reconstructed with negative 
charge sign. For upward going particles, these requirements are reversed. Thus the 
charge sum of all three tracks must be ±1 e. 

• With higher energies, the track momentum resolution and the signal over back- 
ground ratio deteriorate. Thus the total reconstructed momentum must not ex- 
ceed 50 GeV/c. 
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Fig. 6.8: Average energy deposition in the TOF scintillators (top) and average energy depo- 
sition in the last two TOF layers in flight direction (bottom) for data and Monte Carlo electron 
and positron candidates. The Monte Carlo distributions are scaled to the data by integral. 



The requirement for increasing cluster distances within the seed triplets along the flight 
path (see § 16.3.2!) largely distinguishes between downward and upward going particles. To 
make sure the flight direction is correctly recognized, timing information from the TOF 
system is used. Figure shows the distributions of flight times tf for downward and 
upward going particles calculated according to 

tf = {h+t2)/2~h, (6.6) 

where the tj denote the time of passage measured in TOF layer i (i is counted from top to 
bottom). Due to high voltage failures in TOF layer 4 [129] . its timing information is not 
used. The sign of tj in eq. (16.6!) depends on the flight direction. If it is not compatible 
with the flight direction derived from the cluster topology, the event is rejected. The 
time needed by a particle with /5 ^ 1 to pass the distance between the two TOF stations 
is about 4.3 ns. In rare cases particles are attributed flight times considerably smaller 
than this. The reason for these unphysical results is wrong reconstruction of the time of 
passage in the scintillators. To account for the uncertainty of the time measurements, 
is required to be larger than 3.5 ns. 

The energy loss E^cp in a TOF scintillator can be calculated with the Bethe-Bloch 
equation [164] . For an electron or positron in the momentum range 1 to 50GeV/c it 
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amounts to about 2.5 MeV, with fluctuations of about 0.7 MeV |165] . Tlius, to make sure 
that there are three particles traversing the detector, consistent with the signature of a 
converted bremsstrahlung photon, a minimum average energy loss is required in the last 
two TOF layers in the direction of flight. The distribution of this average is shown in 
Figure ES] (bottom). To account for the fluctuations of the energy loss, events with an 
average energy loss smaller than 5 MeV in the last two TOF layers in the direction of 
flight are rejected. 

Apart from protons, nuclei such as He or N have been observed to induce background 
events through hadronic interactions in a few cases. Such particles with Z > 1 deposit 
significantly more energy in the subdetectors than singly charged particles {E^cp > 8 MeV 
for ^He in 1cm scintillator). Thus they can be identified through the average energy 
deposition in all TOF scintillators. The energy loss of an individual particle, due to its 
stochastic nature, follows a Landau distribution [32], which exhibits a long tail towards 
higher values of E^cp- To account for this effect, the truncated mean of the energy 
deposition in the scintillators is calculated by omitting the highest value. Figure 16.81 
(top) shows the distribution of this truncated mean. Events with a mean greater than 
10 MeV are rejected. 
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Fig. 6. 9 : Average of the three highest tracker hit amplitudes for data and Monte Carlo electron 
and positron candidates. The Monte Carlo distributions are scaled to the data by integral. 



In addition to the TOF system, the energy deposition in the silicon strip detectors 
can be used to discriminate against particles with charge g > 1. The amplitudes of the 
hits in the tracker are given in ADC counts, thus the corresponding cut is determined 
from the electron data sample. In Figure 16.91 the distribution of the average of the three 
highest tracker hit amplitudes, regardless of the track to which they belong, is shown. 
The accumulation of positively charged events at higher amplitude is caused mostly by 
helium nuclei. Thus the average amplitude is required not to exceed 110 ADC counts. In 
case of negative total charge, the events at higher amplitude mostly result from wrongly 
reconstructed low momentum background protons. 
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Figf. 6.10: The distribution of the invariant mass at the conversion vertex for data (dots) and 
Monte Carlo (shaded histograms): data events with negative total charge compared to electron 
Monte Carlo (top) and data events with positive total charge compared to proton Monte Carlo 
(bottom). The Monte Carlo distributions are scaled to the data by integral. 



6.4.2 Suppression of Background 

For the suppression of background, use is made of the fact that bremsstrahlung and photon 
conversion imply small opening angles of the particles at the vertices (see § 16.2.11) . while 
in hadronic reactions the particles are emitted under large angles (see § 16.2.21) . However, 
the latter are defined in the center-of-mass frame but observed in the laboratory frame, 
whose relative speed varies with the particles' incident momentum. Thus, in order to 
make the angles independent of the frame of reference, the corresponding invariant mass 
is calculated according to 

= 2 ■ El ■ ^2 ■ (1 - cos , (6.7) 

where 6, Ei and E2 denote the opening angle and the energies of the primary particle and 
the photon and of the conversion pair, respectively. 
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Fig. 6.11: The distribution of the invariant mass at the bremsstrahlung vertex for data 
(dots) and electron and proton Monte Carlo (shaded histograms): data events with negative 
total charge compared to electron Monte Carlo (top) and data events with positive total charge 
compared to proton Monte Carlo (bottom). The Monte Carlo distributions are scaled to the 
data by integral. 



The distributions of invariant mass at the conversion and bremsstrahlung vertex are 
shown for data and electron and proton Monte Carlo in Figures [6. 101 and 16. Hi For events 
with negative total charge, which represent a largely clean electron sample, they reveal 
a narrow shape with a peak at zero. This is in agreement with Monte Carlo results. 
In case of events with positive total charge, consisting of positrons and background, the 
distributions also show a peak at zero, but have a long tail towards higher invariant masses 
which is largely caused by the proton background. The excess due to the positrons is 
clearly visible near zero invariant mass. 

Figure 16.121 shows the distributions for data and for Monte Carlo proton and electron 
events in the invariant mass plane. In agreement with Monte Carlo, data events with 
negative total charge (electrons) accumulate around zero invariant masses. Caused by the 
positrons, such an accumulation is also observed in the distribution of positively charged 
events, where the uniformly distributed proton background indicated by the corresponding 



80 



6.4 Analysis 



Data {negat i ve eve nts} ^ Mont e Car lo electrons 



0.02 ' 



bo 




60 
SO 

40 I 
30 B 

- 10 



w 



0.D2 0.04 Q.Oe 0.03 -D.l 

.31 



-■0 




^ I ^ ' ' I 



140 
120 
ICQ 

40 

■0 



m 



Jnv.conv 



0,02 0.04 Q,04 Q.l 



m 



ln.v.ecinv' 



Data {positive events) 



Monte Carlo protons 




0.02 0.04 0.06 0.08 0. 

..I 



m 



0.0 J O.W 0.Q6 0.03 0,1 

[GeV/c=] 



Ifllf.LCrfl*' 



i^zg. 6.12: Distributions in the invariant mass plane for data events with negative {top left) 
and positive {bottom left) total charge and for Monte Carlo electron {top right) and proton 
{bottom right) events. 



Monte Carlo events is also seen. In order to discriminate against background events, 
restrictions are applied on the invariant masses. These cuts are parameterized as ellipses 
in the invariant mass plane, centered at zero and with half axes in units of the standard 
deviation of the respective distribution in the electron data sample. Events outside the 
ellipses are rejected. In order to keep the positron selection efficiency high, both cut values 
have been set to 2 a. 

6.4.3 Effect of the Mir Space Station 

As stated in § 15. 3^ parts of the Mir space station were within the AMS-01 field of view 
during the four day Mir docking phase of the fiight. Cosmic rays interacting with the Mir 
generated secondary particles |145] . some of which may also have traversed the AMS-01 
detector. Not being of cosmic origin, these secondaries must be rejected in the analysis. 
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positron candidates from this analysis; also indicated is the solid angle obstructed by the Mir 
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Figure 16.131 (left ) displays the incoming directions of downward going particles with 
Z = —1, projected on an x-y-plane above the detector, which were recorded during the 
Mir docking period |145] . 6 and denote the polar and azimuthal angles of the incoming 
particles. There is an apparent excess of particles in a clearly defined area on the left 
side of the plot. It has been shown [145] that these excess particles originated in the 
Mir and were mostly spallation pions and muons from their decay. Hence, since electrons 
and positrons can also result from pion decay, the Mir must be regarded as a possible 
background source. 

In Figure 16.131 (right), the same display is shown for downward going electron and 
positron candidates from this analysis, together with the solid angle obstructed by the 
Mir. No significant excess of particles is observed in the affected solid angle. However, 
particles in the Mir solid angle that were recorded during the docking phase were not 
considered in this analysis. 

6.4.4 Geomagnetic Cutoff 

Energy spectra of cosmic rays are modulated by the geomagnetic field (see § 14.31) . De- 
pending on the direction of incidence and the geomagnetic coordinates of the entry point 
into the magnetosphere, particles with momenta below a certain cutoff are deflected by 
the magnetic field and cannot reach the Earth's proximity. Hence, below geomagnetic 
cutoff the particles detected by AMS-01 must originate from within the magnetosphere. 
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6.5 Correction for Irreducible Background 

They were mostly produced as secondaries through hadronic interactions and trapped 
inside the Earth's radiation belts. 

To discriminate against these secondaries, particle trajectories were individually traced 
back from their measured point of incidence, angle and momentum through the geo- 
magnetic field by numerical integration of the equation of motion using an enhanced 
Runge-Kutta iteration procedure [166] . For this purpose, a reference model of the Earth's 
external field [167] was used in combination with the IGRF-6 coefficients [168j of the 
internal field according to eq. fl4.6l) . 




Fig. 6.14 '- Reconstructed paths of positron candidates in the Earth's magnetic field. The left 
display shows a 16.8GeV/c primary particle, the right one a 2.4GeV/c trapped particle which 
is rejected as being of secondary origin. 

Figure 16.141 illustrates the fiight paths of particles with charge +le by means of two 
examples. The left display shows the reconstructed path of a particle with a momentum of 
16.8 GeV/c, as it approached the Earth from outside the magnetosphere and was detected 
by the AMS-01 experiment over central Northern America. Apparently, this particle is 
of extraterrestrial origin and thus accepted for further analysis. On the right side of 
Figure [U.141 the path of a trapped secondary particle with 2.4 GeV/c momentum is shown 
which exhibits the typical superposition of a helix trajectory with hemisphere bouncing 
and westward drift. It is certain that this particle stems from an interaction with the 
atmosphere near one of the numerous turning points of its trajectory, thus the early part 
of the displayed path is purely hypothetical. Consequently, the particle is not considered 
in the analysis. In general, particles are rejected as secondaries if their trajectory touched 
the surface of the Earth at least once. Particles always staying within a distance of 25 
Earth radii or not crossing the magnetopause are considered as trapped and are also 
rejected. 

6.5 Correction for Irreducible Background 

As seen in Figure I6.12[ the distribution of protons in the invariant mass space does not 
vanish in the signal region around zero invariant mass. The same applies to the back- 
ground from wrongly reconstructed electrons. Consequently, a small fraction of back- 
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ground events will not be rejected by the cut on the invariant masses. This remaining 
irreducible background has to be corrected for, which is accomplished by using Monte 
Carlo data. 

In principle the approach is to run the analysis on an adequate number of proton and 
electron Monte Carlo events as if they were data, determine the amount and momentum 
distribution of particles that are misidentified as positrons, and subtract these from the 
raw positron counts obtained from data. However, such an implicit comparison of Monte 
Carlo and data requires the adjustment of several properties of the simulated events, 
which have not been affected by the geomagnetic field and whose input spectrum is not 
exactly equal to the true fluxes. 

As discussed in § 16. 4^ the geomagnetic field shields the Earth's vicinity from low energy 
particles. However, the geomagnetic cutoff cannot be calculated individually for Monte 
Carlo particles, since their four vector is not defined with respect to the geomagnetic 
coordinates. To correct for the shielding effect, the livetime function T{p), described in 
§ 17.21 and shown in Figure I7.6[ is used. The livetime function gives the effective measure- 
ment time as a function of momentum for singly charged particles. After normalization 
to a maximum value of 1, resulting in the function T(p), its value at a given momentum 
denotes the probability for a singly charged particle to penetrate the geomagnetic field. 
Hence, it serves as a weight for distributions of any event variable in Monte Carlo, parti- 
cularly for the momentum distribution of Monte Carlo background events. The livetime 
function must be evaluated using the reconstructed momentum of the incident particle to 
treat Monte Carlo and real data on an equal footing. 

The incident momentum spectrum of the Monte Carlo particles follows a distribution 
(pMcip) = and therefore differs significantly from the true spectra which have spectral 
indices between 2.7 and 3.4 and are affected by solar modulation. Therefore, the Monte 
Carlo event variables have to be reweighted, since they are correlated with the incident 
momentum. Using the parameterized flux 4>d{p) of protons measured by AMS-01 [16], 
the spectral reweighting function w{p) is calculated as w{p) = (Pd{p)/4>mc{p)- In contrast 
to T{p), w{p) must be evaluated using the incident particle's simulated momentum and 
is not normalized. 

The livetime function as well as the spectral reweighting function are then applied as 
weights in histogramming an event variable x from the proton Monte Carlo sample, mean- 
ing that w{ps) ■ T{pt-) ■ X is histogrammed instead of x, with the simulated momentum ps 
and the reconstructed momentum p^. In particular, the functions correct for the shape of 
the momentum distribution of background protons calculated from Monte Carlo. Subse- 
quently, since w{p) is not normalized, the background distribution must be scaled to the 
data to obtain the total number of background events to be subtracted from the measured 
number of positrons. 

Figure 16.151 illustrates the scaling of the proton Monte Carlo to the positively charged 
data events using the sidebands of the invariant mass distributions. The sidebands are 
defined as the ranges of invariant mass above certain thresholds in which the positron 
contribution to the data is negligible. The thresholds are determined from the electron 
distribution, which is identical to that of positrons. They are set to 0.16GeV/c^ for the 
conversion vertex and to 0.2 GeV/c^ for the bremsstrahlung vertex. The excess in the data 
due to the positron contribution is apparent below threshold. Two scaling factors, one 
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Fig. 6.15: Invariant mass distributions at the conversion (top) and bremsstrahlung vertex 
(bottom) for positively charged events from data and proton Monte Carlo. The Monte Carlo 
distributions have been scaled to the data using the sidebands whose borders are indicated by 
the vertical lines. 
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Fig. 6.16: Distribution of the reconstructed momentum of electrons from data and Monte 
Carlo. The latter have been corrected for the livetime T(p) but not for the distribution of 
momentum. 



from each of the distributions at the two vertices, are determined as the sideband ratios 
of the number of events from data and Monte Carlo. The average of these two factors, 
which differ less than 8 % from each other, is then applied to the reweighted Monte Carlo 
proton background distribution, resulting in the final proton background correction. 

In contrast to protons, the electron data sample is clean, without contamination from 
other particles. Therefore, it can be used directly to reweight the spectrum of the electron 
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Monte Carlo events in order to estimate the background in the positron data resulting 
from wrongly reconstructed electrons. Figure [B.16I shows the distributions Z?e-d(p) and 
-Dg- mc(p) of the measured momenta of electrons from data and Monte Carlo respectively, 
where the latter have been weighted with the normalized livetime function T{p). In a first 
step, the electrons from Monte Carlo that have been wrongly reconstructed as positrons 
are grouped in nine bins according to their simulated momentum following the binning 
in Figure 16.161 For each group, the distribution of measured momenta is separately 
histogrammed. Finally, each of these histograms is scaled with the ratio -Dg- d/-De- mc in 
the respective bins. The sum of the nine resulting histograms then gives the amount and 
distribution of electron background events to be subtracted from the positron sample. 

In total, 1026 electrons and 119 positrons have been found prior to the background cor- 
rection in this analysis. Figure [6. 171 shows the momentum distribution of the uncorrected 
positron sample and the total background correction as a function of momentum, itemized 
into contributions from protons and wrongly reconstructed electrons, which amount to 
24.9 and 6.4 events, respectively. 
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Fig. 6.17: Momentum distribution of positron candidates including background, and the 
background contribution from protons and wrongly reconstructed electrons. 
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7 Results of the Positron Measurement 

7.1 Positron Fraction 

The positron fraction j [e^ + e~) is calculated from the electron and corrected positron 
counts for each momentum bin. It is shown in Figure mi and compared with earlier results 
and with a model calculation based on purely secondary positron production. Table 17.11 
summarizes the results of this analysis. The contributions to the error on the positron 
fraction are discussed in the following. 



momentum 
GeV/c 


range 
GeV/c 








positron 
fraction 


Cstat 




'^sys,a 


1.25 


1-1.5 


11 


3 


3 


0.21 


+0.11 
-0.1 


±0 


±0 


1.75 


1.5-2 


31 


5 


4.8 


0.133 


+0.064 
-0.051 


+0.002 
-0 


±0.006 


2.5 


2-3 


85 


11 


10.7 


0.112 


+0.034 
-0.031 


+0.001 
-0.003 


±0.004 


3.75 


3-4.5 


186 


18 


15.8 


0.078 


+0.021 
-0.018 


+0.001 
-0.003 


±0.004 


5.25 


4.5-6 


172 


16 


10 


0.055 


+0.025 
-0.022 


+0.006 
-0.007 


±0.001 


7.45 


6-8.9 


198 


18 


9 


0.043 


+0.029 
-0.017 


+0.01 
-0.004 


±0.004 


11.9 


8.9-14.8 


195 


22 


14.5 


0.069 


+0.03 
-0.014 


+0.01 
-0.002 


±0.006 


20.6 


14.8-26.5 


109 


20 


15.4 


0.124 


+0.038 
-0.03 


+0.009 
-0.003 


±0.007 


38.2 


26.5-50 


39 


6 


2.9 


0.07 


+0.075 
-0.034 


+0.01 
-0.01 


±0.007 



Tah. 7. 1 : The positron fraction and its momentum dependence, n^- , ng+ j. and ^ denote 
the number of electrons, raw (uncorrected) and corrected positrons, respectively. The last 
three columns give the statistical errors Ustat and the systematic errors (Tsys,b and Usys^a due to 
background and acceptance correction. 



7.1.1 Statistical Errors 

Due to the complexity of the positron fraction analysis, having two sources of irreducible 
background and low statistics, a Bayesian approach based on Monte Carlo simulation has 
been chosen for the determination of the statistical errors |169l 1170] . The aim is to obtain 
the probability distribution of all possible values of the positron fraction which, together 
with the background, can lead to the observed number of particle counts. Given this 
distribution, the confidence intervals are derived by numerical integration. The Monte 
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Fig. 7.1: The positron fraction e~^/{e~^ + e ) measured in this analysis, compared with 
earher results from AMS-01 [2], the combined results from HEAT-e=^ and HEAT-pbar [7] and 
with a model calculation for purely secondary positron production (dashed line). The 

overall errors are given by the outer error bars, while the inner bars represent the systematic 



contributions Usys to the error. 



Carlo procedure outlined below is carried out separately for each of the momentum bins 
defined in Table 17.11 

• In a first step, two random floating point numbers S'e- true and S'e+,true are gener- 
ated by sampling a uniform distribution. They represent the expectation values for 
the probability distributions of "true" numbers of electrons and positrons without 
background. 

• Subsequently, the number of background events are added to S'e+ true- The corrected 
and scaled background counts Bp and B^- from proton and electron Monte Carlo, as 
displayed in Figure 16.171 are the result of a counting experiment and therefore have 
statistical uncertainties, which transform into systematic errors when added to the 
"true" signal. These errors follow Poisson distributions, whose expectation values 
are, to best knowledge, assumed to be equal to Bp and B^,-. Hence, to account for 
their uncertainty. Bp and B^- are smeared according to the individual distributions 
before adding them to S'e+ true- However, the width of the Poisson distributions is 
determined by the statistics of the raw number of Monte Carlo background events 
considered and not by Bp and B^.- , which are the result of correcting and scaling the 
raw counts to the data. As a consequence, the raw Poisson width is altered by the 
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Fig. 7.2: The normalized simulated distributions of the positron fraction /true (shaded his- 
tograms), the fitted functions Dj (blue solid lines), and the resulting confidence interval limits 
(red dashed lines) for each of the nine momentum bins. 



corresponding scale factor. To account for this effect, the raw background counts are 
first smeared according to their respective Poisson distribution and then multiplied 
with the ratio of corrected to raw counts. The resulting numbers of background 
events are denoted by Bp and B^,- . 

• The background counts thus obtained are added to the true number of positrons, 
resulting in an expectation value T+ = 5*6+ true + Bp + B^- for the measured total 
number of positively charged particles. On the other hand, the electron sample is 
largely background- free and T_ = S'e-,true- 

• Two integer random numbers and are finally generated according to Poisson 
distributions with expectation values and T_, respectively. These represent one 
simulated measurement of the number of positron and electron candidates including 
background. If and N_ are exactly equal to the counts actually observed in 
the experiment, the "true" positron fraction /true = 5'e+ true/ ('S'e+ true + 'S'e- true) is 
accepted for further analysis, and the whole procedure is repeated. 
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Figure 17.21 shows the distributions of the simulated positron fraction ftme for each of 
the nine momentum bins. Normahzed to an integral of 1, they represent the probability 
density functions of the positron fraction measured in this analysis, taking into consider- 
ation the amount of background calculated from Monte Carlo. As apparent in the figure, 
the distributions are in all cases parameterized well by a function Df = G{f) ■ L{f), the 
product of a Gaussian and a Landau distribution function G{f) and L{f), respectively. 
The slight asymmetries, which become manifest in the tails towards higher values of /true 
and make the Landau function L{f) necessary, are due to the statistical uncertainty of 
the background correction. 

By numerical integration, the smallest interval is found over which the integral of the 
function Df equals 0.683, yielding the lower and upper limits of the corresponding la 
Gaussian confidence intervals for the appropriate momentum bin. These limits are stated 
as the statistical errors on the positron fraction and are given in Table 17.11 



7.1.2 Systematic Errors 

Since the positron fraction is a ratio of particle fluxes, most sources of systematic error 
such as detector acceptance or trigger efficiency naturally cancel out. The only sources of 
error that must be considered are those which are asjTiimetric with respect to the particle 
charge. 

The background correction is applied to the sample of positron candidates only and 
is therefore a source of systematic error. To a certain degree, the description of the 
experimental setup may be inaccurately implemented in the Monte Carlo program. Fur- 
thermore, in contrast to the production of charged pions, background processes involving 
neutral pion production imply photoconversion with typically low angles between tracks 
emerging from the vertices. Hence, the distribution of invariant masses depends on the 
angular cross sections of charged and neutral pion production. Possible inaccuracies in 
the implementation of the cross sections in the Monte Carlo program must therefore be 
considered. 

The systematic error (Tsys,b from the background correction can be estimated by eval- 
uating the mean deviation of the scaled Monte Carlo background from the data in the 
invariant mass plane. Figure [7751 shows the number density of positively charged particles 
from data ndata (left) and Monte Carlo protons n^c (middle) in the region of invariant 
masses minv > 0.33GeV/c^ at both vertices where the positron signal is negligible. The 
Monte Carlo proton histogram has been scaled to the data using the ratio J2 n-data/ J2 '^mc 
in the bins outside the signal region. The histogram in the right panel of Figure 17.31 shows 
the distribution of relative differences 6,- of data and Monte Carlo bin contents given by 

Or = W ■ (7.1) 

for 7^ with weights w. Although the binning has been chosen rather coarse to 
minimize statistical fluctuations, the remaining statistical uncertainty is accounted for 
by weighting each entry with the inverse of the relative statistical error w = Ut/cTt with 
w ~ in the Poisson limit and Ut = (ridata + 'T-mc)/2. Finally, for display purposes, 
has been scaled to a maximum value of 1. The RMS deviation of (5r leads to a systematic 
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error estimate of about 20 % of the background events. This value is then propagated to 
the positron fraction for each momentum bin, resulting in the asymmetric values of Csys^b 
as given in Table 17.11 No separate systematic error for the background from wrongly 
reconstructed electrons is calculated, since their contribution to the total background is 
small compared to that from protons. 
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-Fig. 7.3: The distribution of particles in the invariant mass space outside the positron signal 
region for positively charged particles from data (left) and scaled Monte Carlo protons (middle). 
The numbers in the histograms denote the particle counts in the respective bins; the distribution 
of the weighted relative difference of the background in data and Monte Carlo (right). 



As a consequence of the East West Effect |171j . in combination with the asymmetric 
layout of the AMS-01 tracker, the product of the detector acceptance and the livetime as 
functions of the direction of incidence of a particle may vary for positrons and electrons, 
though no deviation of their average livetimes is apparent (see section [772]) . This effect 
is accounted for with a second contribution o"sys,a to the systematic error of the positron 
fraction. It is estimated from the mean variation of the difference in livetime of positrons 
and electrons over the detector acceptance. After propagation to the positron fraction, 
the systematic error due to the East West Effect is well below 10 % for all momentum 
bins, except for the highest momenta above 26 GeV, where it amounts to approximately 
14 % of the positron fraction value. 

The total errors atot displayed in Figure 17.11 are calculated according to a^^^ = a'^^^^ + 
"^sys.b + ^sys,a5 Separately for the upper and lower ranges. Similarly, the systematic contri- 
bution (Jsys to the total error is given by a'^^^ = cjgyg + a'^^^^. Clearly, the accuracy of the 
positron fraction is statistically limited by the small number of particle counts. 

7.1.3 Variation of the Positron Fraction with the Value of the Cut 

Figure 17.41 shows the average positron fraction and its total error o"tot as a function of the 
invariant mass cut parameter (see § 16. 4p . The fraction is calculated as the weighted mean 
over the six momentum bins from 3 to 50 GeV/c, which have been corrected for a consid- 
erable amount of background, with the total errors (Ttot as weights. For values of the cut 
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between approximately 1.5 a and 2.5 a, the mean positron fraction remains largely stable, 
showing that the background correction procedure compensates well for the additional 
background events entering the positron sample when the cut is loosened. Beyond 2.5 a, 
the average positron fraction rises constantly, probably due to systematic uncertainties of 
the background correction. However, the cut of 2 a used in this analysis, indicated by the 
vertical line in Figure 17.41 is well within the region in which the background correction 
procedure works reliably. 






0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5 

cut [ct of electron distribution] 



Fig. 7.4 '- Average value of the positron fraction (large dots) and its total error (small dots 
and shaded area) as functions of the cut parameter. The actual cut is indicated by the vertical 
line. 



7.2 Flux Calculation 

As a crosscheck to the measurement of the positron fraction, the absolute incident fluxes 
of electrons and positrons are calculated. The electron flux is then compared to mea- 
surements by other experiments and to the results obtained previously by AMS-01. To 
calculate the fluxes it is necessary to determine the acceptance as well as the livetime for 
positrons and electrons in the context of this analysis. 

The term livetime means the effective amount of time during which cosmic ray particles 
coming from outer space have the opportunity to reach and trigger the experiment. The 
livetime may be significantly lower than the total operation time of the experiment, a 
fact for which four different effects are largely responsible. First, due to a variety of 
reasons, there were several loss-of-signal periods during the flight with average duration 
of roughly one hundred seconds and these must be subtracted from the total operating 
time. Furthermore, the total time during which the trigger system was busy (see § 15.21) 
must also be considered as dead time. 

Besides these rather simple effects, the influence of the Earth and its magnetic field 
causes the livetime to be reduced even in a momentum dependent way. The body of the 
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Earth will obstruct particles arriving from the "wrong" side, except if their momentum is 
low enough so that their trajectories may be bent around the Earth by the geomagnetic 
field. In addition the geomagnetic field will force the trajectories of incoming particles on 
a helix and capture or shield the Earth from particles with momentum below the geomag- 
netic cutoff. These effects depend on the direction of incidence and on the momentum of 
the particle, as well as on the orientation and position of the Space Shuttle, and therefore 
on time. In principle, this means that the experiment may be effectively shut down for 
the measurement of particles in a certain momentum range at a given time. 

From the particle count N{p, 6, 0) in a particular momentum bin p of width Ap and 
knowing the detector acceptance A{p, 6, 0), and the livetime T(j9, 6, 0), one can calculate 
the differential flux as follows: 

d$(p,g,0) N{p,e,<p) 

dp A{p,9,(P)-T{p,e,4>)-Ap- 

If - as is the case with the flux of downward going particles in AMS-01 - the livetime 
is only weakly dependent on direction, the angular distribution of the particle count will 
follow that of the acceptance. Then, one can approximate eq. (17.21) by 

d<l>{p) ^ N{p) 
dp A{p) ■ T{p) - Ap' ^ ■ ^ 

The determination of the detector acceptance and the calculation of the livetime is de- 
scribed in the following two sections. 

7.2.1 Detector Acceptance 

The detector acceptance for the bremsstrahlung conversion process is calculated from 
Monte Carlo, separately for downward and upward going electrons and positrons. In the 
Monte Carlo simulation, particles are emitted above or below the detector from a square 
surface S with a side length of 3.9 m. With Ut the total number of particles emitted from 
5* into the hemisphere facing the detector with an isotropic angular distribution, and Uc 
the number of reconstructed events remaining after all cuts, the acceptance as a function 
of incident momentum is |172j 

A{p) = S.7..^. (7.4) 

Figure [73] shows A{p) separately for upward or downward going electrons and positrons. 
It is of the order of a few cm^- sr and reaches a fiat maximum at approximately 20 GeV/c. 
The decline towards higher momenta is caused by the detector resolution. According 
to the cross section eq. (16.11) . high momentum incident positrons or electrons radiate 
bremsstrahlung photons with higher energy, resulting in higher momentum secondary 
particles which undergo a smaller deflection in the magnetic field of the AMS-01 tracker. 
Therefore, their trajectories are closer to each other and the decreasing cluster separation 
approaches the resolution limit of the silicon strip detectors. By contrast, at low mo- 
mentum secondary particles may be deflected so that they generate multiple separated 
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Fig. 7.5: The acceptance A{p) for downward and upward going positrons and electrons, when 
identified through bremsstrahlung conversion. 



hits in single TOF planes. In this case an event is rejected by the level 3 trigger logic 
of the experiment. Furthermore, for secondary particles with too low momentum, the 
probability that they will not be properly reconstructed increases, leading to a decrease 
of the acceptance. 

Additional material (the ATC, the Space Shuttle's payload bay floor and the support 
structure of AMS-01) is traversed by upward going particles before they enter the detector 
(see § 15.1.61) . Consequently, the acceptance for upward going particles is generally higher 
than for downward going ones. Except for the ATC, the amount of this additional material 
can only be estimated, since precise information is not available. Independent of their 
flight direction, no significant difference in the acceptance for electrons and positrons is 
observed. 



7.2.2 Calculation of Livetime 

For the calculation of the lifetime T as a function of momentum and of the direction of the 
incoming particles a method has been developed which is derived from the determination 
of the geomagnetic cutoff as described in sec. 16.4.41 The basic principle is to determine 
the cutoff rigidity in short time steps for the whole duration of the Space Shuttle flight 
and to count the livetime only for the momentum range above the cutoff. 

In a first step, the acceptance region of AMS-01 is divided into nine bins of equal 
size in cos 6' in the interval [0.7, 1] or [—1, —0.7] for downward or upward going particles, 
respectively, and into eight bins in 0. Furthermore, the momentum range between 1 GeV/c 
and 50GeV/c is divided into eight bins. Then, for every four seconds during the flight 
of the Discovery, her recorded position and attitude are obtained and for each of the 576 
(p, dfi)-bins, a virtual charged particle is started with the corresponding values inside 
the aperture of the detector and propagated backward through the geomagnetic field. If 
the virtual particle fulfills the criteria of a primary cosmic ray particle as described in 



94 



7.2 Flux Calculation 



x10 



I 450 
■§ 400 
350 
300 
250 
200 
150 
100 
50 




Average livetime 

• Downward positrons 
O Downward electrons 
■ Upward positrons 
□ Upward electrons 




3 4 5 6 7 8 910 



20 30 40 

p [GeV/c] 



Fig. 7.6: The livetime, averaged over the detector acceptance, as a function of momentum for 
downward and upward going positrons and electrons. 



section 16.41 the four-second time interval during which the particle has been detected 
is added to the total livetime of the corresponding momentum and acceptance bin after 
subtraction of the trigger dead time. This approach is as accurate as it is CPU-intensive, 
with a total of approximately 5000 GHz-CPU days spent for the calculations. 

The livetime, averaged over the detector acceptance, is displayed in Figure 17.61 for 
downward and upward going positrons and electrons. As a consequence of the rigidity 
cutoff, it is a monotonically increasing function of momentum and converges to approxi- 
mately 140 (55) hours for downward (upward) going particles. These values correspond 
to 75 % (30 %) of the total data taking time of about 184 hours. Since the Space Shuttle 
was mostly oriented towards the zenith (see section 1^751) . and due to obstruction by the 
Earth's body, the livetime for downward going particles exceeds that of upward going 
ones. Concerning the average livetime, no significant difference between positively and 
negatively charged particles is apparent. 



7.2.3 Positron and Electron Fluxes 

Precise information about the thickness of the Space Shuttle's payload bay floor is not 
accessible. Hence, the amount of material underneath the detector, which strongly affects 
the acceptance for upward going particles, cannot be ascertained to a degree which is 
suitable for flux measurements, but only be estimated. As a consequence, in this analysis 
particle fluxes are calculated solely for downward going particles. 

In Figure 17771 the fluxes of downward going positrons and electrons, calculated according 
to eq. (17. 3p . together with results published earlier by AMS-01 [2], and HEAT-e^ [3], are 
displayed together with their statistical errors. The fluxes are in very good agreement with 
previous measurements over the full momentum range, except for a slight discrepancy in 
the electron fluxes between 2 and 3 GeV. Here, at low momenta in combination with low 
statistics, additional inaccuracies of the backtracing through the geomagnetic field are 
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expected to become the dominant source of systematic error to the fluxes. However, this 
effect cancels out in the ratio giving the positron fraction. 
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Over the past decades it has become evident that luminous matter amounts only to a 
small fraction of the energy density in the universe. More than 75 % is accounted for by 
what is called the dark energy, and about 20 % must exist in the form of some kind of 
non-relativistic dark matter. At the same time new physics beyond the Standard Model of 
particle physics is expected to appear when present experimental limitations are exceeded. 
Among the many hypotheses anticipating these new phenomena, supersymmetry yields 
the most promising dark matter candidate, the neutralino. Thus cosmology, astrophysics 
and particle physics may simultaneously approach the same problem and their synergy 
could finally lead to a fundamentally new understanding of basic physics. 

The general purpose detectors, CMS and ATLAS, at the LHC will be able to detect 
sparticles over large ranges of the supersymmetric parameters. Measurements of cosmic 
ray particles, especially of antiparticles such as positrons, can impose strong constraints 
on the nature of new physics beyond the Standard Model. However, cosmic ray positron 
measurements are experimentally very challenging due to the vast proton background. 
While new cosmic ray experiments are under construction or have recently started taking 
data, existing data from past experiments can be re-analyzed using new analysis tech- 
niques. This thesis describes a novel approach of positron identification applied to the data 
of the AMS-01 experiment, namely through the detection of bremsstrahlung conversion 
in a silicon microstrip detector. In contrast to earlier single-track analyses, this approach 
involves the selection and reconstruction of multi-track events. As the discussion of the 
signal process shows, bremsstrahlung from protons is suppressed by a factor of more than 
3 ■ 10^ with respect to positrons due to the dependence of the cross section on the mass of 
the radiating particle. The background to the positron sample, which dominantly stems 
from protons undergoing hadronic interactions and from wrongly reconstructed electrons, 
can largely be suppressed using the topological and geometrical properties of the events. 

In order to obtain the highest positron selection efficiency possible, novel combinato- 
rial track finding algorithms were developed, particularly optimized for the signature of 
converted bremsstrahlung. These algorithms require no additional highly granular subde- 
tector for the determination of track seeds and their momentum resolution is comparable 
to that of single-track events for momenta above 10 GeV. By applying restrictions on the 
invariant mass of particles the background to the positron sample can largely be elimi- 
nated. The remaining background contamination was determined from large samples of 
Monte Carlo data including the effects of the geomagnetic field. It amounted to 26 % of 
the positron counts and was corrected for. No indication was found for additional back- 
ground produced in interactions with the Mir space station. Finally, in order to remove 
atmospheric secondaries from the positron and electron samples, the precise method of 
trajectory backtracing was applied individually to all positron and electron candidates. 

The results of the positron measurement are presented in chapter seven. It is shown 
that the bremsstrahlung approach extends the sensitivity range of AMS-01 to positron 
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momenta up to 50GeV/c, which is far beyond the original scope of the experiment. 
The precision of the positron fraction measurement is clearly statistically limited by the 
small number of particle counts. In the momentum range 1 to 8 GeV the fraction is in 
good agreement with model predictions for background from purely secondary positron 
production, while at higher momenta there is indication for a positron overabundance as 
already reported by other experiments [7j. However, the results of the present analysis are 
subject to different systematic errors, since they were obtained with a different detector 
and a new method of analysis. Furthermore, in contrast to the earlier balloon-borne 
experiments, AMS-01 was operated in space well above the Earth's atmosphere. Thus 
the AMS-01 data lend further weight to the hints of a positron overabundance seen in [7]. 

In addition to the positron fraction, the absolute fluxes of positrons and electrons 
were calculated from the event samples of this analysis. For this purpose, a method was 
developed which allows the determination of the livetime as a function of momentum and 
direction of incidence with high accuracy. The results of the flux calculation are found 
to be consistent with earlier data [21 [3]. In particular, the electron flux for momenta 
above 6 GeV is well in agreement with that obtained from measurements of single-track 
events in the AMS-01 data, confirming the good performance of particle selection with 
the bremsstrahlung approach. 

The positron fraction results from this analysis have been combined with earlier re- 
sets from AMS-01 [2j, HEAT-e± and HEAT-pbar [7], CAPRICE [5] and TS93 [1] (see 
Figure [HH]). For this, a procedure for combining data with asymmetric errors was em- 
ployed [173] . First, all data were grouped in 13 momentum bins. For each data point i 
the corresponding likelihood function Ci was approximated by a Gaussian function with 
a variable standard deviation at. The Ci were parameterized using the values of the 
measurements and their upper and lower errors. The combined likelihood function C, 
from which the combined results and errors are finally determined, is then defined by the 
product of the Ci for a particular momentum bin. The central momentum value for each 
bin is given by the mean momenta of the data in the individual bins. 

The combined positron fraction is shown in Figure 18.21 (top) and is compared to a 
model calculation for the background expected from purely secondary positron production 
(without diffuse reacceleration) [lllj . Table [HTT] summarizes the results. The data indicate 
a positron overabundance for momenta above 6 GeV. The bottom display of Figure 18.21 
shows the cumulative significance, i.e. the significance of all data for momenta p < p' 
with respect to the background-only prediction. The significance was calculated from 
the corresponding probability of the P^r degree of freedom. Because of the possible 
effect of solar modulation on the positron fraction at low momenta, the data point below 
0.3 GeV was not considered in the calculation. Apart from statistical fluctuations on 
the 2 standard deviation level, the positron fraction is in good agreement with purely 
secondary positron production for momenta below 6 GeV. However, when the data at 
higher momenta are considered, the significance of the positron overabundance increases 
to 5.3 standard deviations, which would be reduced to 4.2 standard deviations without 
the results of the present analysis. Therefore, a statistical fluctuation causing the positron 
overabundance in the data with respect to the background-only prediction can now be 
excluded. 
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i^ig. ^.i: The data used for the calculation of the combined positron fraction from: this 
analysis, AMS-01 (single tracks) [2], HEAT-e± and HEAT-pbar [7], CAPRICE [5] and TS93 0]. 
The data were grouped within intervals which are denoted by vertical lines. The solid curve 
shows a model prediction for purely secondary positron production [lllj . 



As stated in § 14.41 besides neutralino annihilation other possible contributions to the 
cosmic ray particle fluxes have been proposed, such as hypothetical Kaluza-Klein particle 
states or galactic pulsars. In order to establish through the analysis of cosmic ray spectra 
if annihilations of dark matter neutralinos are in fact responsible for the observed gamma 
and antiparticle overabundances and to narrow down the parameters of supersymmetry, 
more precise measurements remain necessary. Especially in the case of the cosmic ray 
positron fraction, measurements up to about 300 GeV will allow to precisely distinguish 
possible dark matter positron sources from other contributions. 

In June 2006 the PAMELA satellite experiment [174] . which is designed to study 
charged particles in cosmic rays with a particular focus on antiparticles, was launched 
for an at least three year long mission. PAMELA features, amongst other subdetectors, 
a silicon microstrip tracker in the 0.43 T field of a permanent magnet, a calorimeter and 
a neutron detector. However, with its small geometrical acceptance of 21.5cm^sr, it is 
unlikely that PAMELA will be able to deliver precise data at highest energies relevant to 
a possible dark matter signal. The successor to the AMS-01 experiment, AMS-02 [175] . 
is currently under construction and scheduled for operation on the International Space 
Station (ISS) from 2009 on. Developed on the basis of the AMS-01 design, AMS-02 will 
feature strongly enhanced particle identification by means of a calorimeter, a ring imaging 
Cerenkov detector and a transmission radiation detector (TRD), as well as an improved 
silicon microstrip tracker in the 0.87 T field of a superconducting magnet. With a geo- 
metrical acceptance of the tracking system of 4500 cm^ sr and 900 cm^ sr for the tracker in 
combination with the TRD and the calorimeter, exceeding that of PAMELA by a factor 
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Fig. 8.2: The combined positron fraction of the results from this analysis and earlier results. 
The solid curve shows a model prediction for purely secondary positron production (top); 
the cumulative significance (the significance of all data for momenta p < p') of the positron 
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of more than 40, AMS-02 will conduct cosmic ray spectroscopy with unprecedentedly high 
precision. 
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Tab. 8.1: The combined positron fraction and the data used for its calculation from: this work, AMS-01 (single tracks) [2], HEAT-e^ 
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